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^ ; ABSTRACT 

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and 
Qh' stellar mass as a function of concentration index C r and morphological type in the 

q SDSS. (Our size estimate, a simple analog of the SDSS cmodel magnitude, is new: it 

■ is computed using a combination of seeing-corrected quantities in the SDSS database, 
' and is in substantially better agreement with results from more detailed bulge/disk 
, decompositions.) We also quantify how estimates of the fraction of 'early' or 'late' 

type galaxies depend on whether the samples were cut in color, concentration or light 
profile shape, and compare with similar estimates based on morphology. Our fits show 
that ellipticals account for about 20% of the r-band luminosity density, pL r , and 25% 
• of the stellar mass density, p* ; including SOs and Sas increases these numbers to 33% 

Qs, \ and 40%, and 50% and 60%, respectively. The values of pL r and p*, and the mean 

f^*) . sizes, of E, E+S0 and E+SO+Sa samples are within 10% of those in the Hyde & 

Bernardi (2009), C r > 2.86 and C r > 2.6 samples, respectively. Summed over all 
galaxy types, we find p* ~ 3 x lO 8 Af Mpc~ 3 at z ~ 0. This is in good agreement with 
expectations based on integrating the star formation history. However, compared to 
' most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 

c"^ , at M* ~ 5 x 10 Mq. The stellar mass density further increases at large masses if we 

assume different IMFs for elliptical and spiral galaxies, as suggested by some recent 
chemical evolution models, and results in a better agreement with the dynamical mass 
function. 

■ We also show that the trend for ellipticity to decrease with luminosity is primarily 

because the E/S0 ratio increases at large L. However, the most massive galaxies, A/* > 
5x 10 11 M Q , are less concentrated and not as round as expected if one extrapolates from 
lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation 
histories with recent radial mergers. Finally, we show that the age-size relation is flat 
for ellipticals of fixed dynamical mass, but, at fixed Afdyn, SOs and Sas with large 
sizes tend to be younger. Hence, samples selected on the basis of color or C r will yield 
different scalings. Explaining this difference between E and SO formation is a new 
challenge for models of early-type galaxy formation. 

Key words: galaxies: formation - galaxies: haloes - dark matter - large scale structure 
of the universe 



1 INTRODUCTION 

Each galaxy has its own peculiarities. Nevertheless, even 
to the untrained eye, sufficiently well-resolved galaxies can 
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be separated into three morphological types: disky spirals, 
bulgy ellipticals, and others which are neither. The mor- 
phological classification of galaxies is a field that is nearly 
one hundred years old, and sample sizes of a few thousand 
morphologically classified galaxies are now available (e.g. 
Fukugita et al. 2007; Lintott et al. 2008). However, such 
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eyeball classifications are prohibitively expensive in the era 
of large scale sky surveys, which image upwards of a few 
million galaxies. Moreover, the morphological classification 
of even relatively low redshift objects from ground-based 
data is difficult. Thus, a number of groups have devised au- 
tomated algorithms for discerning morphologies from such 
data (e.g. Ball et al. 2004 and references therein). 

In parallel, it has been recognized that relatively simple 
criteria, using either crude measures of the light profile (e.g. 
Strateva et al. 2001), the colors (e.g. Baldry et al. 2004), or 
some combination of photometric and spectroscopic infor- 
mation (Bernardi et al. 2003; Bernardi & Hyde 2009) allow 
one to separate early-type galaxies from the rest. Because 
they are so simple, these tend to be more widely used. The 
main goal of this paper is to show how samples based on 
such crude cuts compare with those which are based on 
the eyeball morphological classifications of Fukugita et al. 
(2007). We do so by comparing the luminosity, stellar mass, 
size and velocity dispersion distributions for cuts based on 
photometric parameters with those based on morphology. 
These were chosen because the luminosity function is stan- 
dard, although it is becoming increasingly common to com- 
pare models with <f>(M*) rather than cj>(L) (e.g. Cole et al. 
2001; Bell et al. 2003; Panter et al. 2007; Li & White 2009); 
the size distribution 4>(R) has also begun to receive con- 
siderable attention recently (Shankar et al. 2009b); and the 
distribution of velocity dispersions 4>(a) (Sheth et al. 2003) is 
useful, amongst other things, to reconstruct the mass distri- 
bution of super-massive black holes (e.g. Shankar et al. 2004; 
Tundo et al. 2007; Bernardi et al. 2007; Shankar, Weinberg 
& Miralda-Escude 2009; Shankar et al. 2009a) and in studies 
of gravitational lensing (Mitchell et al. 2005). 

Section [2] describes the dataset, the photometric and 
spectroscopic parameters derived from it, and the subsam- 
ple defined by Fukugita et al. (2007). This section shows 
how we use quantities output from the SDSS database to 
define seeing-corrected half-light radii which closely approx- 
imate the result of bulge + disk decompositions. We describe 
our stellar mass estimator in this section as well; a detailed 
comparison of it with stellar mass estimates computed by 
three different groups is presented in Appendix [X] The re- 
sult of classifying objects into two classes, on the basis of 
color, concentration index, or morphology are compared in 
Section [3] Luminosity, stellar mass, size and velocity dis- 
persion distributions, for the Fukugita et al. morphological 
types are presented in Section |3J where they are compared 
with those based on the other simpler selection cuts. This 
section includes a discussion of the functional form, a gen- 
eralization of the Schechter function, which we use to fit 
our measurements. We find more objects with large stellar 
masses than in previous work (e.g. Cole et al. 2001; Bell 
et al. 2003; Panter et al. 2007; Li & White 2009); this is 
the subject of Section [5] where implications for the match 
with the integrated star formation rate, and the question of 
how the most massive galaxies have evolved since z ~ 2 are 
discussed. 

While we believe these distributions to be interesting in 
their own right, we also study a specific example in which 
correlations between quantities, rather than the distribu- 
tions themselves, depend on morphology. This is the cor- 
relation between the half-light radius of a galaxy and the 
luminosity weighted age of its stellar population. Section [6] 



shows that the morphological dependence of this relation 
means it is sensitive to how the 'early-type' sample was se- 
lected, potentially resolving a discrepancy in the recent lit- 
erature (Shankar & Bernardi 2009; van der Wei et al. 2009; 
Shankar et al. 2009c). A final section summarizes our results, 
many of which are provided in tabular form in Appendix [B] 
Except when stated otherwise, we assume a spatially 
flat background cosmology dominated by a cosmological 
constant, with parameters (f2 m , S1a) = (0.3, 0.7), and a Hub- 
ble constant at the present time of Ho = 70 km s _1 Mpc _1 . 
When we assume a different value for Ho, we write it as 
H = WOh km s _1 Mpc _1 . 



2 THE SDSS DATASET 
2.1 The full sample 

In what follows, we will study the luminosities, sizes, velocity 
dispersions and stellar masses of a magnitude limited sample 
of ~ 250, 000 SDSS galaxies with 14.5 < m Pe trosian < 17.5 
in the r— band, selected from 4681 deg 2 of sky. In this band, 
the absolute magnitude of the Sun is M r ,0 = 4.67. 

The SDSS provides a variety of measures of the light 
profile of a galaxy. Of these the Petrosian magnitudes and 
sizes are the most commonly used, because they do not de- 
pend on fits to models. However, for some of what is to fol- 
low, the Petrosian magnitude is not ideal, since it captures 
a type-dependent fraction of the total light of a galaxy. In 
addition, seeing compromises use of the Petrosian sizes for 
almost all the distant lower luminosity objects, leading to 
systematic biases (see Hyde & Bernardi 2009 for examples). 

Before we discuss the alternatives, we note that there 
is one Petrosian based quantity which will play an impor- 
tant role in what follows. This is the concentration index 
C r , which is the ratio of the scale which contains 90% of 
the Petrosian light in the r-band, to that which contains 
50%. Early-type galaxies, which are more centrally concen- 
trated, are expected to have larger values of C r , and two 
values are in common use: a more conservative C r > 2.86 
(e.g. Nakamura et al. 2003; Shen et al. 2003) and a more 
cavalier C r > 2.6 (e.g. Strateva et al. 2001; Kauffmann et al. 
2003; Bell et al. 2003). We show below that from the first 
approximately two-thirds of the sample comes from E+S0 
types, whereas the second selects a mix in which E+S0+Sa's 
account for about two-thirds of the objects. 

The SDSS also outputs deV or exp magnitudes and 
sizes which result from fitting to a deVaucouleur or expo- 
nential profile, and f racDev, a quantity which takes values 
between and 1, which is a measure of how well the deVau- 
couleur profile actually fit the profile (1 being an excellent 
fit). In addition, it outputs cmodel magnitudes; this is a 
very crude disk+bulge magnitude which has been seeing- 
corrected. Rather than resulting from the best-fitting lin- 
ear combination of an exponential disk and a deVaucouleur 
bulge, the cmodel magnitude comes from separately fitting 
exponential and deVaucouleur profiles to the image, and 
then combining these fits by finding that linear combina- 
tion of them which best-fits the image. Thus, if m cxp and 
mdcv are the magnitudes returned by fitting the two models, 
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Figure 1. Comparison between apparent magnitudes and sizes obtained from performing full bulge+disk decompositions (denoted 
deV+exp), with those output by or constructed from parameters in the SDSS database. In all cases, symbols with error bars show the 
mean relation and the error on the mean, and dashed lines show the actual rms scatter. Filled circles in left panels show results for the 
SDSS cmodel magnitudes and effective radii (see text for details) and right panels are for the SDSS deV (red stars; the effective radius 
is the value in the SDSS database multiplied by ^/b/a, i.e. %/ba) or exp (blue open squares; the effective radius is the value in the SDSS 
database a) quantities. Cyan triangles in bottom left panel show the result of picking either the deV or exp size, based on which of the 
corresponding magnitudes were closer to the cmodel magnitude. Although these triangles are almost indistinguishable from the filled 
circles, the rms scatter is substantially larger, particularly at small f racdev. 



then 

10 -0.4m adss _ cmodel = ^ _ fracdeV ) 10 -0.4m Elp 

+ f racdeV 10"°' 4mdoV . (1) 

Later in this paper, we will be interested in seeing- 
corrected half-light radii. We use the cmodel fits to define 
these sizes by finding that scale r eiCmod ei where 

iQ-0.4m Bdss _ Model /"^sdss — cmodel 

— = (1 - fracdeV) 2vr / d8 6 I exp (9) 

2 Jo 

f' r sdss — cmodel 

+ f racdeV 2tt / dO 9 Jdev(0),(2) 

Jo 

where I is the surface brightness associated with the two fits. 
Note that the SDSS actually performs a two dimensional fit 
to the image, and it outputs the half-light radius of the long 
axis of the image a, and the axis ratio b/a. The expression 
above assumes one dimensional profiles, so we use the half- 



light radius a of the exponential fit, and y/b/a x a = %/ba 
for the deVaucouleur fit. We describe some tests of these 
cmodel quantities shortly. 

We would also like to study the velocity dispersions of 
these objects. One of the important differences between the 
SDSS-DR6 and previous releases is that the low velocity dis- 
persions (a < 150 km s _1 ) were biased high; this has been 
corrected in the DR6 release (see DR6 documentation, or 
discussion in Bernardi 2007). The SDSS-DR6 only reports 
velocity dispersions if the S/N in the spectrum in the rest- 
frame 4000 — 5700 A is larger than 10 or with the status 
flag equal to 4 (i.e. this tends to exclude galaxies with emis- 
sion lines). To avoid introducing a bias from these cuts, we 
have also estimated velocity dispersions for all the remain- 
ing objects (see Hyde & Bernardi 2009 for more discussion). 
These velocity dispersions are based on spectra measured 
through a fiber of radius 1.5 arcsec; they are then corrected 
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Figure 2. Top panels: Comparison of cmodel magnitudes (right) and sizes (left) with those obtained from performing full bulge+disk 
decompositions as function of cmodel sizes. The sample was divided in three bins based on the shape of the light profile (as indicated in 
the panels). Thin solid curves (blue, green and red) show fits to equations (0 and Q. Except for the sample with fracDev > 0.8, the 
coefficients of these fits are given in Table [T] For fracDev > 0.8, the coefficients in Table [T] are based on the larger sample of Hyde & 
Bernardi (2009) (see text for the origin of the small offsets); this results in the thick (magenta) solid curve shown. Bottom panels: Similar 
to panels on left of Figure[T] but with cmodel magnitudes and sizes corrected following equations J3]l and J4}. In all panels, symbols show 
the mean relation, error bars show the error on the mean and dashed lines (bottom) show the rms scatter. 



to r e /8, as is standard practice. (This is a small correction.) 
The velocity dispersion estimate for emission line galaxies 
can be compromised by rotation. In addition, the disper- 
sion limit of the SDSS spectrograph is 69 km s - , so at 
small a the estimated velocity dispersion may both noisy 
and biased. We will see later that this affects the velocity 
dispersion function. The size and velocity dispersion can be 
combined to estimate a dynamical mass; we do this by set- 
ting M dyn = 5R e a 2 /G. 

2.2 A morphologically selected subsample 

Recently, Fukugita et al. (2007) have provided morphologi- 
cal classifications (Hubble type T) for a subset of 2253 SDSS 
galaxies brighter than m Pe t = 16 in the r— band, selected 
from 230 deg 2 of sky. Of these, 1866 have spectroscopic in- 
formation. Since our goal is to compare these morphological 
selected subsamples with those selected based on relatively 



simple criteria (e.g. concentration index), we group galaxies 
classified with half-integer T into the smaller adjoining in- 
teger bin (except for the E class; see also Huang & Gu 2009 
and Oohama et al. 2009). In the following, we set E (T — 
and 0.5), SO (T = 1), Sa (T = 1.5 and 2), Sb (T = 2.5 and 3), 
and Scd (T = 3.5, 4, 4.5, 5, and 5.5). This gives a fractional 
morphological mix of (E, SO, Sa, Sb, Scd) = (0.269, 0.235, 
0.177, 0.19, 0.098). Note that this is the mix in a magnitude 
limited catalog - meaning that brighter galaxies (typically 
earlier-types) are over-represented. 

2.3 cmodel magnitudes and sizes 

As a check of our cmodel sizes, we have performed deVau- 
couleurs bulge + exponential disk fits to light profiles of a 
subset of the objects; see Hyde & Bernardi (2009) for a de- 
tailed description and tests of the fitting procedure. If we 
view these as the correct answer, then the top left panel 
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Table 1. Coefficients used in equations {3} and J1J to correct 
sizes and magnitudes for sky subtraction problems. 
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0.3 < fracDev < 0.8 
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0.001 
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of Figure [T] shows that the cmodel magnitudes are in good 
agreement with those from the full bulge+disk fit, except 
at fracDev~ and fracDev~ 1 where cmodel is fainter by 
0.05 mags (top left). This is precisely the regime where the 
agreement should have been best. As discussed shortly, the 
discrepancy arises mainly because the SDSS reductions suf- 
fer from sky subtraction problems (see, e.g., SDSS DR7 doc- 
umentation), whereas our bulge-disk fits do not (see Hyde 
& Bernardi 2009 for details). Comparison with the top right 
panel shows that cmodel is a significant improvement on 
either the deV or the exp magnitudes alone. 

The bottom panels show a similar comparison of the 
sizes. At intermediate values of fracDev, neither the pure 
deVaucouleur nor the pure exponential fits are a good de- 
scription of the light profile, so the sizes are also biased 
(bottom right). However, at fracDev=l, where the deVau- 
couleurs model should be a good fit, the deV sizes returned 
by the SDSS are about 0.07 dex smaller than those from 
the bulge+disk decomposition. There is a similar discrep- 
ancy of about 0.05 dex with the SDSS Exponential sizes at 
fracDev=0. We argue below that these offsets are related 
to those in the magnitudes, and are primarily due to sky 
subtraction problems. 

The filled circles in the bottom left panel show that our 
cmodel sizes (from equation [2} are in substantially better 
agreement with those from the bulge+disk decomposition 
over the entire range of fracDev, with a typical scatter of 
about 0.05 dex (inner set of dashed lines). For comparison, 
the triangles show the result of picking either the deVau- 
couleur or exponential size, based on which of these mag- 
nitudes were closer to the cmodel magnitude (this is essen- 
tially the scale that the SDSS uses to compute model colors). 
Note that while this too removes most of the bias (except at 
small fracDev), it is a substantially noisier estimate of the 
true size (outer set of dashed lines). This suggests that our 
cmodel sizes, which are seeing corrected, represent a signif- 
icant improvement on what has been used in the past. 

The SDSS reductions are known to suffer from sky sub- 
traction problems which are most dramatic for large objects 



or objects in crowded fields (see DR7 documentation). The 
top panels in Figure [2] show this explicitly: while there is lit- 
tle effect at small size, the SDSS underestimates the bright- 
nesses and sizes when the half-light radius is larger than 
about 5 arcsec. Note that this is actually a small fraction 
of the objects: 6% of the objects have sdss-cmodel sizes 
larger than 5 arcsec; 13% are larger than 4 arcsec. Whereas 
previous work has concentrated on mean trends for the full 
sample, Figure shows that, in fact, the difference depends 
on the type of light profile - galaxies with fracDev > 0.8 
(i.e. close to deVaucouleur profiles) are more sensitive to 
sky-subtraction problems than later-type galaxies. Some of 
this is due to the fact that such galaxies tend to populate 
more crowded fields. 

To correct for this effect, we have fit low order poly- 
nomials to the trends; the solid curves in the top panels 
of Figure [2] show these fits. Except for the sample with 
fracDev > 0.8, we use these fits to define our final corrected 
cmodel sizes by: 

Te, cmodel — ^e.sdss — cmodel + Cr + Cn r e , sdss — cmodel 

+ Cr 2 r 

e, sdss — cmodel (3) 

and 

'^cmodel — '^sdss — cmodel + Cmo + Cmi r e , sdss-cmodel 

+ Cmi r Ej sdss -cmodel , (4) 

where the coefficients Cmo, Cmi, Cm?, Cro, Cr\ and CY2 
are reported in Table [1] 

For objects with fracDev > 0.8, the trends we see are 
similar to those shown in Fig. 5 of Hyde & Bernardi (2009) , 
which were based on a (larger) sample of about 6000 early- 
type galaxies. The thick solid (magenta) line in the top pan- 
els of Figure [2] show the Hyde-Bernardi trends, with a small 
offset to account for the fact that they did not integrate the 
fitted profiles to infinity (because the SDSS, to which they 
were comparing, does not), whereas we do. The thick solid 
curve differs from the thin one at sizes larger than about 
5 arcsec. Since the thick curve is based on a larger sample, 
we use it to define our final corrected cmodel sizes. The cor- 
rections are again described by equations ((3| and Q, with 
coefficients that are reported in Table [1] 

The bottom panels of Figure [5] show that these cor- 
rected quantities agree quite well with those from the full 
bulge+disk fit, even at small and large fracDev. 

2.4 Stellar Masses 

Stellar masses M» are typically obtained by estimating 
M*/L (in solar units), and then multiplying by the rest- 
frame L (which is not evolution corrected). In the following 
we compare three different estimates of M*. All these esti- 
mates depend on the assumed IMF: Table [2] shows how we 
transform between different IMFs. 

The first comes from Bell et al. (2003), who report 
that, at 2 = 0, log 10 (M*/L r )o = 1.097 (g — r)o + zp, 
where the zero-point zp depends on the IMF (see their Ap- 
pendix 2 and Table 7). Their standard diet-Salpeter IMF 
has zp = —0.306, which they state has 70% smaller M*/L r 
at a given color than a Salpeter IMF. In turn, a Salpeter 
IMF has 0.25 dex more M*/L r at a given color than the 
Chabrier (2003) IMF used by the other two groups whose 
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Table 2. M*/L IMF offsets used in this work. Offset is with 
respect to the Salpeter (1955) IMF: logi M*/L (IMF Salpeter) 
= logio Mt/L (IMF Reference) + Offset 



IMF 


Offset [dcx] 


Refecence 


Kennicut 


0.30 


Kennicut (1983) 


Scalo 


0.25 


Scalo (1986) 


diet-Salpeter 


0.15 


Bell & de Jong (2001) 


pseudo-Kroupa 


0.20 


Kroupa (2001) 


Kroupa 


0.30 


Kroupa (2002) 


Chabrier 


0.25 


Chabrier (2003) 


Baldry & Glazbrook 


0.305 


Baldry & Glazbrook (2003) 



mass estimates we use (see Table[3]for conversion of different 
IMFs). Therefore, we set zp = -0.306+0.15-0.25 = -0.406, 
making 



log 10 (M»Bcii/Lr)o = 1-097 (g - r) - 0.406. 



(5) 



We then obtain M*Beii by multiplying by the SDSS r— band 
luminosity. When comparing with previous work, we usu- 
ally use Petrosian magnitudes, although our final results 
are based on the cmodel magnitudes which we believe are 
superior. 

Note that this expression requires luminosities and col- 
ors that have been k- and evolution-corrected to z = (E. 
Bell, private communication). Unfortunately, these correc- 
tions are not available on an object-by-object basis. Bell 
et al. (2003) report that the absolute magnitudes brighten 
as 1.3z and g — r color becomes bluer as 0.3z. Although 
these estimates differ slightly from independent measure- 
ments of evolution by Bernardi et al. (2003) and Blanton et 
al. (2003), and more significantly from more recent deter- 
minations (Roche et al. 2009a), we use them, because they 
are the ones from which equation ([5]) was derived. Thus, in 
terms of restframe quantities, 



logio 



Mr. 



1.097 ( 5 -r)-0.406-0.4(M r -4.67)-0.19z. 

(6) 



If we use the restframe r — i color and L r luminosity instead, 
then 



logio 



Mr 



-- 1.431 (r-i)-0.122-0.4(M r -4.67)-0.23z. 

(7) 

These two estimates of M» will differ because there is 
scatter in the (g — r) — {r — i) color plane. Unfortunately, Bell 
et al. do not provide a prescription which combines differ- 
ent colors. Although we could perform a straight average of 
these two estimates, this is less than ideal because the value 
of r — i at fixed g — r may provide additional information 
about M*/L. In practice, we will use the g — r estimate as 
our standard, and r — i to illustrate and quantify intrinsic 
uncertainties with the current approach. 

The second estimate of M» is from Gallazzi et al. (2005). 
This is based on a likelihood analysis of the spectra, assumes 
the Chabrier (2003) IMF, and returns M,/L z . The stellar 
mass M^Gaiiazzi is then computed using SDSS petrosian z- 
band restframe magnitudes (i.e., they were fc-corrected, but 
no evolution correction was applied). In this respect, they 
differ in philosophy from M».b c ii, in that the M* estimate 
is not corrected to z = 0. In practice, since we are mainly 



interested in small lookback times from z = 0, for which 
the expected mass loss to the IGM is almost negligible, this 
almost certainly makes little difference for the most massive 
galaxies. These estimates are only available for 205, 510 of 
the objects in our sample (~ 82%). The objects for which 
Gallazzi et al. do not provide stellar mass estimates are lower 
luminosity, typically lower mass objects; we show this explic- 
itly in Figure [201 

A final estimate of M,/L r comes from Blanton & 
Roweis (2007), and is based on fitting the observed colors 
in all the SDSS bands to templates of a variety of star for- 
mation histories and metallicities, assuming the Chabrier 
(2003) IMF. In the following we use the Blanton & Roweis 
stellar masses computed by applying the SDSS petrosian 
and model (restframe) r-magnitudes to these mass-to-light 
estimates: M*p etr o and M, Model - Blanton & Roweis also pro- 
vide mass estimates from a template which is designed to 
match luminous red galaxies; we call these M.lrg- Note 
that in this case M*/L is converted to M* using (restframe) 
model magnitudes only, since the Petrosian magnitude is 
well-known to underestimate (by about 0.1 mags) the mag- 
nitudes of LRG-like objects (Blanton et al. 2001, DR7 doc- 
umentation). To appreciate how different the LRG template 
is from the others, note that it allows ages of upto 10 Gyrs, 
whereas that for the others, the age is more like 7 Gyrs. 

A detailed comparison of the different mass estimates is 
presented in AppendixfS] This shows that to use the Blanton 
& Roweis (2007) masses, one must devise an algorithm for 
choosing between M, Model and M.lrg- In principle, some 
of the results to follow allow one to do this, but exploring 
this further is beyond the scope of the present paper. On 
the other hand, to use the Gallazzi et al. (2005) estimates, 
one must be wary of aperture effects. Finally, stellar masses 
based on the k + e corrected r — i color show stronger system- 
atics than do the g — r based estimates of M* . Therefore, in 
what follows, our prefered mass estimate will be that based 
on k + e corrected cmodel r— magnitudes and g — r colors 
(i.e., equation[6]with restframe and evolution corrected mag- 
nitudes and colors). Note that we believe the cmodel mag- 
nitudes to be far superior to the Petrosian ones. Of course, 
when we compare our results with previous work which used 
Petrosian magnitudes (Section [5}, we do so too. 



3 MORPHOLOGY AND SAMPLE SELECTION 

This section compares a number of ways in which early-type 
samples have been defined in the recent literature, with the 
morphological classifications of Fukugita et al., and discusses 
what this implies for the 'red' fraction. When we show lumi- 
nosities, they have been corrected for evolution by assuming 
that the magnitudes brighten with redshift as 1.3z. 



3.1 Simple measures of the light profile 

Concentration index, axis ratio, and f racdev have all been 
used as proxies for selecting red, massive, early-type galax- 
ies. So it is interesting to see how these quantities correlate 
with morphological type. 

The bottom right panel of Figure [3] shows the distribu- 
tion of all objects in the Fukugita et al. sample in the space 
of concentration versus luminosity. The two horizontal lines 
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Figure 3. Distribution of morphological types from Fukugita et al. (2007) in the space of CV vs luminosity. Bottom right panel shows the 
full sample. Horizontal lines (same in all panels) show two popular cuts for selecting early-types. Text in each panel shows the fraction of 
objects in the panel which lie above these lines. Cyan symbols with error bars and flanked by dashed curves (same in each panel) show 
the median relation and the rms scatter defined by a sample selected following Hyde & Bernardi (2009). 



show C r = 2.6 and CV = 2.86, the two most popular choices 
for selecting early-type samples. The symbols with error bars 
show the mean concentration index at each L if the sample 
is selected following Hyde & Bernardi (2009): i.e. fracDev 
= 1 in g- and r-, and r-band b/a > 0.6. To this we add 
the condition logio(r e , 9 /r ei7 .) < 0.15, which is essentially a 
cut on color gradient (Roche et al. 2009b). This removes a 
small fraction (< 2%) of late- type galaxies which survive 



the other cuts. Dashed curves show the scatter around the 
Cr-luminosity relation in the Hyde & Bernardi sample. 

The other panels show the result of separating out the 
various morphological types. Whereas Es and SOs occupy 
approximately the same region in this space, most Es (93%) 
have C r > 2.86, whereas the distribution of C r for SOs is 
somewhat less peaked. (Here, the percentages we quote are 
per morphological type, in the Fukugita et. al. sample - 
meaning a sample that is magnitude-limited to m ri p e t < 16, 
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Figure 4. Same as previous figure, but as a function of stellar mass rather than luminosity. 



with no 1/Knax weighting applied.) Samples restricted to 
C r > 2.6 have a substantial contribution from both Sas 
(74% of which satisfy this cut) and SOs (for which this frac- 
tion is 93%), and this remains true even if C r > 2.86 (50% 
of Sas and 77% of SOs) . Thus, it is difficult to select a sam- 
ple of Es on the basis of concentration index alone. On the 
other hand, the larger mean concentration of the Hyde & 
Bernardi selection cuts suggest that they produce a sample 
that is dominated by ellipticals/SOs, and less contaminated 
by Sas (see Section [3.51 and Table below) . Figure 2] shows 



that replacing luminosity with stellar mass leads to similar 
conclusions. 

Before moving on, notice that although the mean con- 
centration increases with luminosity and stellar mass in the 
Hyde & Bernardi sample, this is no longer the case at the 
highest L r or M*: we will have more to say about this 
shortly. 

Figures [5] and [6] show a similar analysis of the axis ratio 
b/a. The different symbols with error bars (same in all pan- 
els) show samples selected to have C r > 2.6, C r > 2.86, and 
following Hyde & Bernardi. At M r < — 19 or so, the mean 
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Figure 5. Similar to Figure [3] but now in the space of b/a vs luminosity. Magenta open squares and green diamonds, each with error 
bars and flanked by dashed curves (same in all panels), show the samples selected using cuts in C r larger of 2.6 and 2.86, respectively. 



b/a in the first two cases increases with luminosity upto 
M T — —22.7 or so; it decreases for the brightest objects. At 
low L, the sample with CV > 2.6 has smaller values of b/a 
on average, though the scatter around the mean is large. 
Morevoer, while there are essentially no Es with b/a < 0.6 
about 26% of SOs have b/a < 0.6. On the other hand, a little 
less than half the Sa's and many Scd's also have b/a > 0.6 
(because they are face on). Evidently, just as CV alone is not 
a good way to select a pure sample of ellipticals, selecting 
on b/a alone is not good either. 

The filled red circles (with error bars) in the top left 



panel show that, for Es, b/a « 0.85 independent of M r , ex- 
cept at the highest luminosities where b/a decreases. This 
independence of M r differs markedly from that in either of 
the C r samples, but is reproduced by the Hyde & Bernardi 
sample, for which b/a = 0.8 except at M r < —22.7 where it 
decreases. The difference of about 0.05 in b/a arises because 
the Hyde & Bernardi sample includes some SOs (we quantify 
this in Table [3] below) , for which b/a ~ 0.7 (filled red circles 
in top right panel). This leads to an important point. While 
it has long been known that b/a tends to increase with lumi- 
nosity, even in 'early-type' samples, our Figure [S] shows that 
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Figure 6. Same as previous figure, but as a function of stellar mass rather than luminosity. 



this increase is driven by the changing morphological mix - 
the change from SOs to Es - at M r > —23. Whether this is 
due to environmental or pure secular evolution effects is an 
open question. 

On the other hand, there is a plausible, environmentally 
driven model for the decrease in b/a at the highest L. This 
decrease has been expected for some time (see Gonzalez- 
Garcia & van Albada 2005; Boylan-Kolchin et al. 2006 and 
references therein) - it was first found by Bernardi et al. 
(2008). This is thought to indicate an increasing incidence 
of radial mergers, since these would tend to result in more 



prolate objects. The decrease in concentrations at these high 
luminosities (Figure [3} is consistent with this picture, as is 
the fact that most of these high luminosity objects are found 
in clusters. All of the preceding discussion remains true if one 
replaces luminosity with stellar mass. 

Before concluding this section, we note that we have 
also considered the quantity fracdev which plays an 
important role in the selection cuts used by Hyde & 
Bernardi (2009). The vast majority of ellipticals (85%) have 
fracdev=l, with only a percent or so having fracdev< 0.8. 
The distribution of fracdev has a broader peak for SOs, but 
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Figure 7. Color-magnitude relation for the different morphological types in the Fukugita et al. sample. Solid and dashed curves (same 
in each panel) show the mean of the 'red' and 'blue' sequences, and their thickness, which result from performing double-Gaussian fits 
to the color distribution at fixed magnitude, of the full SDSS galaxy sample (from Bernardi et al. 20f0, in preparation). Dotted green 
line shows the luminosity dependent threshold used to separate 'red' from 'blue' galaxies (equation \Efy . 



they otherwise cover the same range as Es: only 37% have 
fracdev< 1. However, 70% of Sas have fracdev< 1, whereas 
for Scs's, only 10% have f racdev> 0.4. (Note that the above 
percentages were computed in the magnitude limited cata- 
log, i.e. were not weighted by 1/Vm a x-) 



3.2 Colors 

In addition to simple measures based on the light profile, 
or more commonly, as an alternative to such methods, color 
is sometimes used as a way to selecting early types. This is 
typically done by noting that the color-magnitude distribu- 
tion is bimodal (e.g. Baldry et al. 2005), and then adopting 
a crude approximation to this bimodality (e.g. Zehavi et 
al. 2005; Blanton & Berlind 2007; Skibba & Sheth 2009). 
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Figure 8. Same as previous figure, but now for stellar mass in place of luminosity. 



Figure [7] shows this bimodal distribution in the Fukugita 
et al. sample (Figure [8] shows the corresponding color- M, 
relation). The dotted green line in Figure [7] shows 

g - r = 0.63 - 0.03 (M r + 20). (8) 

It runs approximately parallel to the 'red' sequence, and is 
similar to that obtained by subtracting —0.17 mags from 
equation (4) in Skibba & Sheth (2009); it is shallower than 
equation (7) in Skibba & Sheth (2009) or equation (1) of 
Young et al. (2009) (note that we fc-correct to z — and we 
use h = 0.7). 'Red' galaxies are those which lie above this 
line; 'blue' lie below it. 



However, note that many late types (Sb and later) lie 
above this line - these tend to be edge-on disks. In addi- 
tion, some Es lie below it. (See Huang & Gu 2009 for a 
more detailed analysis of such objects, which show either a 
star forming, AGN or post-starburst spectrum.) We intend 
to present a more detailed study of the morphological de- 
pendence of the color-magnitude, stellar mass and velocity 
dispersion relations in a future paper (Bernardi et al. 2009, 
in preparation). For our purposes here, we simply note that 
Figures [7] and [5] illustrate that cuts in color are not a good 
way to select early-type galaxies. 
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Figure 9. Joint distribution of age and stellar mass in the Fukugita et al. (2007) sample. Cyan filled circles, green diamonds and 
magenta squares show the median age at fixed stellar mass for a subsample selected following Hyde & Bernardi (2009), a subsample with 
C r > 2.86, and a subsample with C r > 2.6, respectively. The dashed lines show the lcr range around the median. 



3.3 Ages 

Later in this paper, we will also study correlations between 
stellar age and galaxy mass, size and morphology. The age 
estimates we use, from Gallazzi et al. (2005), are based on a 
detailed analysis of spectral features. Since this is the same 
analysis that provided M*Gaiiazzi, errors in age and stellar 
mass are correlated (see Bernardi 2009 for a detailed dis- 
cussion) : erroneously large M, will tend to have erroneously 
large age as well. 



Figure[5]shows the age-M* correlation for the objects in 
the Fukugita et al. sample for which age estimates are avail- 
able. This shows that, for any given morphological type, 
massive galaxies tend to be older (this correlation is not due 
to correlated errors). However, as expected, the later types 
tend to be substantially younger: Whereas two-thirds of the 
Es in this sample are older than 8 Gyrs, only half of SOs, 
one-quarter of Sas, and fewer than 10 % of the later types 
(Sb, Sc, etc.) are this old. The bottom right panel suggests 
that the age-M, distribution separates into two populations 
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Figure 10. Joint distribution of age and stellar mass in the full sample. Small dots show a random subsample of the galaxies selected 
following Hyde & Bernardi (2009). Cyan filled circles, green diamonds and magenta squares show the median age at fixed stellar mass 
for a subsample selected following Hyde & Bernardi (2009), a subsample with CV > 2.86, and a subsample with C r > 2.6, respectively. 
The dashed lines show the la range around the median. The fraction of galaxies with luminosity weighted ages older than 8 Gyrs and 
younger than 6 Gyrs, for each of the selection methods, are shown. 



- one which is younger than about 7 Gyrs and another which 
is older. However, this is not simply correlated with morpho- 
logical type: the top panels shows that this bimodality is also 
present in the E and SOs samples. 

To study this further Figure [lOl shows the age-M* dis- 
tribution in a random subsample of the galaxies selected fol- 
lowing Hyde & Bernardi (2009) from the full SDSS catalog. 
Note that although 90% of the objects are older than 6 Gyrs, 
this selection clearly includes a population of younger ob- 
jects. This population of 'rejuvenated' early-type galaxies 
has been the subject of some recent interest (e.g. Huang 
& Gu 2009; Thomas et al. 2010). Cyan filled circles show 
that the median age increases with stellar mass and dashed 
lines show the la range around the median. Green diamonds 
and magenta squares show the median age at fixed stellar 
mass for objects with C r > 2.86 and C r > 2.6, respectively. 
At large M, , both these samples produce similar age-M* 
relations to the Hyde-Bernardi sample; at smaller M, , al- 
lowing smaller C r includes younger galaxies. These median 
relations are superimposed on the panels of Figure [9] 



3.4 The red and blue fractions 

There is considerable interest in the 'build-up' of the red 
sequence, and the possibility that some of the objects in the 
blue cloud were 'transformed' into redder objects. We now 
compare estimates of the red or blue fraction that are based 
on color and concentration, with ones based on morphology. 
In particular, we show how these fractions vary as a function 
of L, M, , and a. 



All the results which follow are based on samples which 
are Petrosian magnitude limited, so in all the statistics we 
present, each object is weighted by V^ X (L), the inverse of 
the maximum volume to which it could have been seen. This 
magnitude limit is fainter for the full sample (m Pet = 17.75) 
than it is for the Fukugita et al. subsample (mp et = 16), and 
note that Vmax depends on our model for how the luminosi- 
ties evolve (absolute magnitudes brighten as 1.3z). 

The panels on the left of Figure [11] show how the mix 
of objects changes as a function of luminosity, stellar mass, 
and velocity dispersion, for the crude but popular hard cuts 
in concentration and color (as described in the the previous 
sections). 

Figure [11] shows that the fraction of objects which sat- 
isfies the criteria used by Hyde & Bernardi (2009) increases 
with increasing L, M* and <r, except at the largest values 
(about which, more later). Requiring C r > 2.86 instead re- 
sults in approximately 5% to 10% more objects (compared 
to the Hyde & Bernardi cuts) at each L or M, ; although, in 
the case of a ~ 300 km s _ , this cut allows about 20% more 
objects. Relaxing the cut to CV > 2.6 allows an additional 
15%, with slightly more at intermediate L and M*. 

Selecting objects redder than a luminosity dependent 
threshold (equation [H]) which runs parallel to the 'red' se- 
quence allows even more objects into the sample, but com- 
bining the color cut with one on b/a reduces the sample to 
one which resembles C r < 2.86 rather well. The cut in b/a 
is easy to understand, since edge-on discs will lie redward 
of the color cut even though they are not early-types - the 
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Figure 11. Left panels: The 'red' fraction in the full sample (m rj p e t < 17.5), where being red means the object lies redward of a 
luminosity dependent threshold which runs parallel to the 'red' sequence (orange stars); or has concentration greater than 2.6 (magenta 
squares); or greater than 2.86 (green diamonds); or satisfied the Hyde & Bernardi selection cuts (red crosses). Also shown is the result 
of combining the color cut with one on b/a (brown triangles). Right panels: The fraction of objects of a given L, M*, a and R e , as later 
and later morphological types are added to the Fukugita et al. sample (m r Pet < 16). In all cases, each object was weighted by V^ X (L), 
the maximum volume to which it could have been seen, given the apparent magnitude limit. 



additional cut on b/a is an easy (but rarely used!) way to 
remove them. 

It is interesting to compare these panels with their coun- 
terparts on the right of Figure 111! in which later and later 
morphological types are added to the Fukugita et al. sub- 
sample which initially only contains Es. This suggests that 
the Hyde & Bernardi selection will be dominated by E, 
C r > 2.86 will be dominated by E+SOs, and CV > 2.6 will 
be dominated by E+SO+Sas. We quantify this in the next 
subsection. Figure [12] shows a similar comparison with the 



blue fraction. Note that the contamination of the red frac- 
tion by edge-on discs is a large effect: 60% of the objects 
at log 10 (M*/Mq) = 10.5 are classified as being red, when 
E+SO+Sas sum to only 40%. Figure [13] shows this more di- 
rectly: the reddest objects at intermediate luminosities are 
late-, not early-type galaxies. 

Before concluding this section, we note that both con- 
centration cuts greatly underpredict the red fraction of the 
most luminous or massive objects, as does the application 
of a b/a cut to the straight color selection or the Hyde & 
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Figure 12. Left panels: Same as left hand panel of previous figure, but now showing the 'blue' fraction - the objects which did not qualify 
as being 'red'. Right panels: Similar to right hand panel of previous figure, but now starting with later types and adding succesively 
more early-types. 



Bernardi selection. The most luminous or massive elliptical 
galaxies in the Fukugita et al. sample show the same be- 
havior. I.e., the most massive objects are less concentrated 
for their luminosities than one might have expected by ex- 
trapolation from lower luminosities and their light profile 
is not well represented by a pure deVaucoleur law. This is 
consistent with results in the previous section where, at the 
highest luminosities, b/a tends to decrease with luminosity 
(Figure [5}. These trends suggest an increasing incidence of 
recent radial mergers for the most luminous and massive 
galaxies. 



3.5 Distribution of morphological types in 
differently selected samples 

Much of the previous analysis suggests that the Hyde & 
Bernardi selection will produce a sample that is dominated 
by Es, C r > 2.86 will include more SOs and Sas, and C r > 2.6 
will include Sas and later types. Table [3] shows the distri- 
bution of types in subsamples selected from the Fukugita 
et al. (2007) catalog to have C r > 2.6, C r > 2.86 and fol- 
lowing Hyde & Bernardi (2009). Of the 1596 objects in the 
magnitude limited catalog, 1009, 802 and 470 satisfy these 
cuts. The Table shows that, in samples where C r > 2.6, 
54% of the objects are Sa or later. This fraction falls to 40% 
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Figure 13. Depedence of the bimodality in the color distribution on morphological type, for a few bins in luminosity. Dashed histogram 
shows the distribution of the SDSS sample (m r> p e t < 17.5), while solid histograms show the distribution in the Fukugita et al. sample 
("ir.Pet < 16) as later and later morphological types are added. Note that the reddest objects at intermediate luminosity are late-type 
galaxies. 



Table 3. The morphological mix in differently selected samples, 
from Vma X weighted counts in the Fukugita et al. (2007) sample 
restricted to M r < —19; numbers in brackets are from the raw 
counts. 



Type 


HB09 


C r > 2.86 


Cr > 2.6 


E 


0.69 (0.73) 


0.38 (0.51) 


0.26 (0.43) 


SO 


0.23 (0.20) 


0.22 (0.23) 


0.20 (0.22) 


Sa 


0.07 (0.06) 


0.25 (0.17) 


0.30 (0.20) 


Sb 


0.01 (0.01) 


0.12 (0.07) 


0.19 (0.12) 


Scd 




0.03 (0.02) 


0.05 (0.03) 



for C r > 2.86 and to less than 10% for the Hyde-Bernardi 
cuts (these numbers are obtained after weighting by V^axj 
so they do not depend on the selection effect associated with 
the apparent magnitude limit of the catalog) . 

These numbers indicate that Es comprise at least two- 
thirds of the Hyde-Bernardi sample, but in a sample where 
C r > 2.86, to reach this fraction one must include SOs, and if 
C r > 2.6, then reaching this fraction requires adding Sas as 
well. Stated differently, E's comprise more than two-thirds 
of a Hyde-Bernardi sample, but about one-third of a sample 
with C r > 2.86 and one quarter of a sample with CV > 2.6. 
If we weight each object by its stellar mass, then (E+S0)s 
account for (72+21)% of the total stellar mass in a Hyde & 
Bernardi sample, (47+23)% in a sample with C r > 2.86, and 



(39+23)% if C r > 2.6. These differences will be important 
in Section [6] 



4 DISTRIBUTIONS FOR SAMPLES CUT BY 
MORPHOLOGY OR CONCENTRATION 

We now show how the luminosity, stellar mass, size and 
velocity dispersion distributions - 4>(L), <f>(M*), (j>(R e ) and 
0(cr) - depend on how the sample was denned. We use the 
same popular cuts in concentration as in the previous Sec- 
tions, C r > 2.86 and CV > 2.6, which we suggested might 
be similar to selecting early-type samples which, in addition 
to Es, include SOs + Sas, and SOs + Sas + Sbs, respectively. 
We then make similar measurements in the Fukugita et al. 
subsample, to see if this correspondence is indeed good. 

In this Section, we use cmodel rather than Petrosian 
quantities, for the reasons stated earlier. The only place 
where we continue to use a Petrosian-based quantity is when 
we define subsamples based on concentration, since CV is the 
ratio of two Petrosian-based sizes, or for comparison with re- 
sults from previous work. 

The results which follow are based on samples which 
are Petrosian magnitude limited, so in all the statistics we 
present, each object is weighted by V^^(L), the inverse of 
the maximum volume to which it could have been seen. In 
addition to depending on the magnitude limit (m Pe t < 17.5 
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for the full sample, and m Pe t < 16 for the Fukugita subsam- 
ple), the weight V^ m (L) also depends on our model for how 
the luminosities have evolved. A common test of the accu- 
racy of the evolution model is to see how (1^/Knax) , the ratio 
of the volume to which an object was seen to that which 
it could have been seen, averaged over all objects, differs 
from 0.5. In the full sample, it is 0.506 for our assumption 
that the absolute magnitudes evolve as 1.32; had we used 
1.62z (Blanton et al. 2003), it would have been 0.509. On 
the other hand, if we had ignored evolution entirely, it would 
have been 0.527. 

In addition, SDSS fiber collisions mean that spectra 
were not taken for about 7% of the objects which satisfy 
mpet < 17.5. We account for this by dividing our V^~ ax 
weighted counts by a factor of 0.93. This ignores the fact that 
fiber collisions matter more in crowded fields (such as clus- 
ters); so in principle, this correction factor has some scatter, 
which may depend on morphological type. We show below 
that, when we ignore this scatter, then our analysis of the 
full sample produces results that are in good agreement with 
those of Blanton et al. (2003) , who account for the fact that 
this factor varies spatially. This suggests that the spatial 
dependence is small, so, in what follows, we ignore the fact 
that it (almost certainly) depends on morphological type. 



The peak of Xcf>(X) occurs at X max = X* (a//3) 1/7/3 , where 
C max = ~a.fi. Since a and f3 are usually positive, errors typ- 
ically act to decrease the height of the peak. Since the net 
effect of errors is to broaden the distribution, hence extend- 
ing the tails, errors also tend to decrease f3. The expression 
above shows that, in the O/O, 3> 1 tail, errors matter more 
when j3 is large. 

Fitting to equation (|10fl rather than to equation ((9} is 
a crude but effective way to estimate the intrinsic shape 
(i.e., to remove the effect of measurement errors on the fit- 
ted parameters) , provided the measurement errors are small. 
The rms errors on (In L r , In M» , In R e , In a) are indeed small: 
fcrr = (0.05, 0.25, 0.15, 0.15), and so it is the results of these 
fits which we report in what follows. However, to illustrate 
which distributions are most affected by measurement er- 
ror, we also show results from fitting to equation (|9}; in 
most cases, the differences between the returned parameters 
are small, except when /3 > 1. 

In practice, the fitting was done by minimizing 



£)[ W -]og 10 (]n(I0) 0^(0) 



(12) 



where \ji was log 10 of the V^" ax weighted count in the ith 
logarithmically spaced bin (and recall that, because of fiber 
collisions, the weight is actually 



4.1 Parametric form for the intrinsic distribution 

We will summarize the shapes of the distributions we find 
by using the functional form 



4>(X)AX = </>« 



dX 
~X~' 



(9) 



This is the form used by Sheth et al. (2003) to fit the dis- 
tribution of velocity dispersions; it is a generalization of the 
Schechter function commonly fit to the luminosity function 
(which has ft = 1, a slightly different definition of a). We 
have found that the increased flexibility which /3 7^ 1 al- 
lows is necessary for most of the distributions which follow. 
This is not unexpected: at fixed luminosity, most of the ob- 
servables we study below scatter around a mean value which 
scales as a power-law in luminosity (e.g. Bernardi et al. 2003; 
Hyde & Bernardi 2009). Because this mean does not scale 
linearly with L, and because the scatter around the mean 
can be significant, then if </>(!/) is well-fit by a Schechter 
function, it makes little physical or statistical sense to fit 
the other observables with a Schechter function as well. 



4.2 Effect of measurement errors 

In practice, we will also be interested in the effect of mea- 
surement errors on the shape of the distribution. If the errors 
are Lognormal (Gaussian in lnX) with a small dispersion, 
then the observed distribution is related to the intrinsic one 



by 

Otp(0) 

c 



\nX X<j>{X)p(\nO\\nX) 



04>{0) ( 1 + c 



(a -13 (0/0*fY - f(OIO*f 



(10) 
(11) 



4.3 Covariances between fitted parameters 

When fitting to equation ©, a reasonable understanding of 
the covariance between the fitted parameters can be got by 
asking that all parameter combinations give the same mean 
density (Sheth et al. 2003): 



Px 



!>* X, 



r[(i + a)/fl 

V{fx/P) ■ 



(13) 



In practice, 0* is determined essentially independently of the 
other parameters, so it is the other three parameters which 
are covariant. Hence, it is convenient to define 



{X)=PKI 



(14) 



If j3 is fixed to unity, then this becomes (X) = X* a. But 
if j3 is not fixed, then a further constraint equation can be 
got by asking that all fits return the same peak position or 
height. For distributions with broad peaks, it may be better 
to instead require that the second central moment, 



2 

ox 



(X z ) - (Xy where {X A ) = Xt 



a T[(2 + «)//?] 



Y{a/p) 



(15) 

be well reproduced. Thus, the covariance between a and 
j3 is given by requiring that tix/{X) equal the measured 
value for this ratio. The changes to these correlations are 
sufficiently small if we instead fit to equation (|10l) . so we 
have not presented the algebra here. 



4.4 Distributions for samples cut by concentration 

Figure [14] shows 4>{L) and 0(M») in the full sample (top 
curve in top panel), when one removes objects with CV < 2.6 
(second from top), objects with CV < 2.86 (third from top), 
and when one selects early-types on the basis of a number 
of other criteria (bottom, following Hyde & Bernardi 2009). 
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M r [mag] Log 10 M. [M Sun ] 

Figure 14. The effect of different selection cuts on the cmodel luminosity and stellar mass functions: black crosses, magenta triangles, 
green diamonds and blue squares show the full sample, a subsample with C r > 2.6, a subsample with C r > 2.86, and a subsample 
selected following Hyde & Bernardi (2009). Top panels: Smooth solid curves show the fits to the observed distributions (equation [9{ 
while dashed curves (almost indistinguishable from solid curves) show the intrinsic distributions (equation 1 1 P . Cyan solid line shows 
our fit to the full sample for M r < —20 and M t > lO lo Af0. The parameters of these fits are reported in Tables [BTUB4I Red dotted line 
shows the measurement associated with Petrosian quantities (i.e. from Figures j 191 and 1201 . Bottom panels: Same as top panel, but now 
each set of data points and dashed curve are shown after dividing by their associated solid curve. 



For comparison, the dotted curves show the measurement 
in the full sample when Petrosian quantities are used (i.e. 
from Figures RHJ1 and |2"0"]) . 

The solid curves show the result of fitting to equa- 
tion and the dashed curves (almost indistinguishable 
from the solid ones, except at the most massive end) result 
from fitting to equation ()10[) instead, so as to remove the 
effects of measurement error on our estimate of the shape 
of the intrinsic distribution. To reduce the dynamic range, 
bottom panels show each set of curves divided by the asso- 
ciated solid curve (i.e., by the fit to the observed sample). 
The dotted lines in these bottom panels show that Petrosian 
based counts lie well below those based on cmodel quan- 
tities at M r < -20 or log 10 M*/M > 10.6. The dashed 
lines in the bottom right panel show that the intrinsic dis- 
tribution has been noticably broadened by errors above 
log lo Af»/M > 11. 

Figure [15] shows a similar analysis of <j)(R e ) and </>(er). 
For 4>{cj) we also compare our results with those of Sheth et 
al. (2003). The Sheth et al. analysis was based on a sample 
selected by Bernardi et al. (2003), which was more like CV > 
2.86 at large masses, but because of cuts on emission lines 



and S/N in the spectra, had few low mass objects. Indeed, 
at large a, our measured 4>{o) is similar to theirs. 

Tables [B"T1lB4l rcport the parameters of the fits (to equa- 
tions [5] and I10p shown in Figures [TP] and 1151 It is worth 
noting that, for a given sample, the fits return essentially 
the same value of (j>* in all the tables, even though this was 
not explicitly required. And note that </>(M,) and 4>i (J ) are 
the distributions that are most sensitive to errors; the for- 
mer because the errors are large, and the latter because (5 is 
large. 



4.5 Comparison with morphological selection 

Figures [Pol andfTTl show cj>(L), <f>(M»), 4>{R e ) and 4>{cr) for the 
Fukugita et al. sample as one adds more and more morpho- 
logical types. The smooth curves (same in each panel) show 
the fits to the samples shown in Figures [T4land [T5l In the top 
panels of Figure [16] only, we also show the result of remov- 
ing objects with CV < 2.6 before making the measurements. 
This allows a direct comparison with one of the curves from 
the previous Figure. Notice that this gives results which are 
very similar to those from the larger (fainter, deeper) full 
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Figure 15. Same as previous figure, but now for the size and velocity dispersion. Cyan solid line shows our fit to the full sample for 
1.5 < -Re/kpc < 20 and a > 125 km s _1 . The top right panel also shows the fits obtained by Sheth et al. (2003) to the observed (solid 
grey line) and intrinsic (dashed grey) 4>(cr) distribution. 



sample; the different magnitude limits do not matter very 
much. 

The results of fitting equation (|10|) to the Fukugita sub- 
samples are provided in Tables IB1HB4I which also show the 
associated luminosity and stellar mass densities, and the 
mean sizes. Ellipticals account for about 19% of the r-band 
cmodel luminosity density and 25% of the stellar mass den- 
sity; including SOs increases these numbers to 33% and 41%, 
and adding Sas brings the contributions to 50% and 60% re- 
spectively. 

Note that the number, luminosity and stellar mass den- 
sities of Es - about 10~ 3 Mpc~ 3 , 0.2 x lO 8 L Mpc~ 3 and 
0.6 x lO 8 A/0Mpc -3 respectively - are very similar to those 
of the early-type sample selected following Hyde & Bernardi 
(2009), as is the mean half-light radius of 3.2 kpc. Some of 
this match is fortuitous - we showed before that E's account 
for about 70% of this sample, not 100%. However, this lack 
of purity is balanced by the fact that Hyde & Bernardi se- 
lect about 75% of the Es, not 100%: the purity and com- 
pleteness effects approximately cancel. Similarly, although 
requiring CV > 2.86 produces counts which are similar to 
those of E+S0s, about 40% of the sample is made of Sas 
and later types, but the purity again approximately cancels 
the incompleteness. Finally, the counts when CV > 2.6 are 
similar to those in the E+S0+Sa sample, although 25% of 
the objects are of later type. 



5 THE STELLAR MASS FUNCTION IN THE 
FULL SAMPLE 

Our stellar mass function has considerably more objects at 
large A/* than reported in previous work. Before we quantify 
this, we show the results of a variety of tests we performed 
to check that the discrepancy with the literature is real. This 
is important, since the high mass end has been the subject 
of much recent attention (e.g., in the context of the build-up 
of the red sequence). 



5.1 Consistency checks 

We first checked that we were able to reproduce previous re- 
sults for the luminosity function. These have typically used 
Petrosian rather than cmodel magnitudes, and different Hq 
conventions. The top panel in Figure IT51 shows the result of 
estimating the luminosity function using Petrosian magni- 
tudes (using the Vmax method and code as in the previous 
sections) and the curve show the Schechter function fits re- 
ported by Blanton et al. (2003). This agreement shows that 
our algorithms correctly transform between different Ho con- 
ventions, and between different definitions of fc-corrections. 
It also shows that varying the evolution correction between 
the value reported by Blanton et al. (1.6z) and that from Bell 
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Figure 16. Distributions for objects of different morphological types in the Fukugita et al. (2007) sample (cyan filled circles with error 
bars). A subsample with CV > 2.6 is also shown in the 'E' panels (red triangles). Smooth solid curves are same as in Figurt tl4l and [T5l 
from top to bottom in each panel, they show the full SDSS sample (black), a subsample with CV > 2.6 (magenta), a subsample with 
C r > 2.86 (green), and a subsample selected following Hyde & Bernardi (2009) (blue). 
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Figure 17. 



Continued. 
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Figure 18. Luminosity function in the 01 r— and r— bands (top 
and bottom panels), determined from r— band Petrosian appar- 
ent magnitudes in the range 14.5 — 17.5, for the two choices of 
the pure luminosity evolution parameter advocated by Blanton 
et al. (2003) and Bell et al. (2003): 1.6* and 1.32, respectively. 
Dashed line in top panel shows the fit reported by Blanton et al. 
(2003); solid line in bottom panel shows that reported by Bell et 
al. (2003). 



et al. (1.32), which we use when estimating stellar masses, 
makes little difference. 

The bottom panel shows the Schechter function fit re- 
ported by Bell et al. (2003); there is good agreement. How- 
ever, note that here we have not shifted the Petrosian mag- 
nitudes brightwards by 0.1 mags for galaxies with C r > 2.6 
(Bell et al. did this to account crudely for the fact that Pet- 
rosian magnitudes underestimate the luminosity of early- 
types). At faint luminosities, the measurements oscillate 
around the fits, suggesting that fits to the sum of two 
Schechter functions will provide better agreement, but we 
do not pursue this further here. At the bright end, we find 
slightly more objects than either of the Bell et al. (top) or 
Blanton et al. (bottom) fits, but a glance at Fig. 5 in Blanton 
et al. shows that the fit they report slightly underestimates 
the counts in the high luminosity tail. 

In contrast to the good agreement for the luminosity 
function, our estimates of the stellar mass function (Fig- 
ure I20|) show a significant excess relative to previous work 
at M* > 1O 11,5 M0. (Note that, to compare with previous 



Figure 19. Luminosity function in the r— band determined from 
Petrosian apparent magnitudes in the range 14.5 < m r ,pet < 
17.5. Red line shows the fit reported by Bell et al. (2003). Cyan 
line shows our fit to equation l(9} for log 10 M^/Mq > 10.5, with 
parameters reported in Table [4] 



work, we convert from M t / L r to M„ by using the Petrosian 
magnitude. However, we do not shift the Petrosian magni- 
tudes brightwards by 0.1 mags for galaxies with C r > 2.6 
- a shift that was made by Bell et al. (2003). If we had 
done so the excess would be even larger). To ensure that 
the discrepancy with previous fits is not caused by outliers, 
we removed galaxies with L > 10 11 L Q or M* > 10 115 M Q 
which differ by more than 0.3 dex from the linear fit (solid 
line in Figure I2ip , in luminosity or M* . Except for the g — r 
based M, , where we show results before and after removal 
of these outliers, all the other measurements shown in Fig- 
ure[20]are from the sample in which these outliers have been 
removed. Note that while removing outliers makes the plot 
slightly cleaner (see Figure [20}, the discrepancy at high M 4 
remains. 

The discrepancy is most severe if we use stellar masses 
from Gallazzi et al. (2005). (Note that they do not pro- 
vide stellar mass estimates for fainter, typically lower mass 
objects.) The discrepancy is slightly smaller if we use equa- 
tion ([6]) to translate g — r color into M*/L r (for which both 
g— r and L r are evolution corrected; if we had used restframe 
quantities without correcting for evolution the discrepancy 
would be even worse). Using r — i instead (equation!?} yields 
results which are more similar to the original Bell et al. 
(2003) fit. And finally, (/>(M*) based on Af*p c tro of Blanton 
& Roweis (2007) has the lowest abundances of all. 

The dotted line shows that our measurement of the dis- 
tribution of M*p c t ro is well-fit (except for a small offset) by 
the formula reported by Li & White (2009) , which was based 
on their own estimate of <f)(M*) from Blanton & Roweis 
M*Petro- The fact that we find good agreement with their 
fit suggests that our algorithm for estimating <f)(M t ) from a 
given list of M» values is accurate. 

The dashed-dot-dot green line shows the fit reported 
by Panter et al. (2007) who computed stellar masses for a 
sample of 3 x 10 5 SDSS galaxies based on the analysis of the 
spectral energy distribution of the SDSS spectra. While this 
fit lies slightly below our data at high M t , the discrepancy 
is smaller than it is for most of the other fits we show. 
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Figure 20. Stellar mass functions estimated from Petrosian magnitudes. Top: Green open diamonds show the distribution of M, 
estimated from g — r colors (equation [6} for the full SDSS sample; filled black circles show the same after removing outliers at L > 10 Lq 
or M, > 10 115 Mq (see Figure [21}, and filled grey circles use r — i (equation [7} instead. Grey open squares use M, estimates from 
Gallazzi et al. (2005). Open red triangles show Af„,p ctro from Blanton & Roweis (2007). Solid cyan line shows our fit to equation (01 with 
parameters reported in Table [4] for log 10 M*/M@ > 10.5. Solid grey curve shows the fit reported by Bell et al. (2003), dashed-dot blue 
curve that of Cole et al. (2001), dashed-dot-dot green curve that of Panter et al. (2007), and dotted red curve the fit from Li &; White 
(2009), all transformed to H = 70 km s^Mpc" 1 and Chabrier (2003) IMF. Bottom: Same as top panel, but now all quantities have 
been normalized by the Bell et al. (2003) fit (and the results corresponding to the green open diamonds are not shown). 
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Figure 21. Petrosian L r vs M* (M* computed using equation[6}. 
Small dots show a representative subsample of the galaxies. Solid 
green line shows the linear fit measured for M, > 10 11 Mq. 
Dashed and dotted green lines show the 1 and 2 — a range around 
the median, respectively. Red filled circles show outliers that were 
removed when determining <j>(L) and <f>(M*). 



We argued previously that aperture effects may have 
inflated the Gallazzi et al. masses slightly. As a check, we 
recomputed masses from the g — r color using equation ©, 
but now, using the Fiber color output by the SDSS pipeline. 
In contrast to the model color, which measures the light on a 
scale which is proportional to the half-light radius, this mea- 
sures the light in an aperture which has the same size as the 
SDSS fiber. The spectro-photometry of the survey is suffi- 
ciently accurate that this is a meaningful comparison (e.g. 
Roche et al. 2009a). Note that this gives abundances which 
are larger than those based on the model colors. Moreover, 
they are almost indistinguishable from those of M*Gaiiazzi- 

The cyan solid curve shows the result of fitting equa- 
tion to our measurements of (j>(M t ) based on g — r color 
(we do not show fits to equation 1101 because none of the 
other fits in the literature account for errors in the stellar 
mass estimate) when log 10 M,/M > 10.5. The best-fit pa- 
rameters are reported in Table [3] The Table also reports 
the fit to the full sample (i.e. log 10 M«/M > 8.5). (These 
differ from those reported in the previous section, because 
here they are based on Petrosian rather than cmodel mag- 
nitudes.) The Figures show the results for M T < —20 and 
log 10 M*/Mq > 10.5 because this is the regime of most in- 
terest here. 

Our estimates depart from the Bell et al. fit at densities 
of order 10 -4 Mpc -3 . However, their analysis was based on 
only 412 deg 2 of sky, for which the expected number of ob- 
jects on the mass scale where we begin to see a discrepancy 
(M* ~ 10 Mq) is of order tens. This, we suspect, is the 
origin of the discrepancy between our results and theirs - 
we are extrapolating their fit beyond its regime of validity. 

Bell et al. (2003) did not account for errors, so the most 
straightforward way to quantify the increase we find is to 
compare our measured counts with their fit. If we do this for 
our Petrosian-based counts, then the stellar mass density in 
objects more massive than (1, 2, 3) x 10 11 Mq is ~ (3, 24, 86) 
percent larger than one would infer from the Bell et al. fit. 
However, Bell et al. attempted to account for the fact that 
Petrosian magnitudes underestimate the total light by shift- 



ing galaxies with C r > 2.6 brightwards by 0.1 mags. Since 
we have not performed such a shift, the appropriate com- 
parison with their fit is really to use our measured cmodel 
based counts, so the difference between our counts and Bell 
et al. are actually larger. We do this in the next section. 

If we compare our estimate of the stellar mass den- 
sity in objects more massive than (1,2,3) x 10 11 Mq with 
those from the Li & White (2009) fit, then our values are 
~ (109,202,352) percent larger. 

5.2 Towards greater accuracy at large M* 

It is well known that the cmodel luminosities are more re- 
liable at the large masses where the discrepancy in <f>(M*) 
is largest. Therefore, Figure [22] compares various estimates 
of 4>{M- t ) based on cmodel magnitudes. In this case, the es- 
timates based on M*Modci of Blanton & Roweis (2007) pro- 
duce the lowest abundances, (but note they are larger than 
those based on M*p ctr o in Figure [217)) . whereas those based 
on M*lrg are substantially larger, as one might expect (c.f. 
Figure [X3|. The M*lrg abundances are also in good agree- 
ment (slightly smaller) with those based on g — r color (equa- 
tion [6] with cmodel magnitudes), except at smaller masses. 
(Although it is not apparent because of how we have chosen 
to plot our measurements, above 10 11 Mq, the abundances 
based on M*Gaiiazzi and M*lrg are in good agreement.) 

We argued previously that we believe the g — r masses 
are more reliable than those based on r — i. Therefore, we 
only show fits to the distribution of g — r derived masses: 
solid and dashed curves show fits to equation (0 and (|1U|1 
respectively (the latter account for broadening of the dis- 
tribution due to errors in the determination of M*). Both 
result in larger abundances than the observed abundances 
based on Petrosian quantities - here represented by the fit 
shown in the previous figure - although the intrinsic distri- 
bution we determine for the cmodel based masses is similar 
to the observed distribution of Petrosian-based masses. 

If we sum up the observed counts to estimate the stellar 
mass density (M* from equation[SJ) in objects more massive 
than (1, 2, 3) x 10 11 Mq, then the result is ~ (30, 68, 170) per- 
cent larger than that one infers from the Bell et al. (2003) 
fit. Using our fit to the observed distribution (values be- 
tween round brackets in Table[B2] for log 10 M* /Mq > 10.5) 
gives similar results: stellar mass densities ~ (21,60,160) 
percent larger than those from the Bell et al. (2003) fit. 
In practice, however, the Bell et al. fit tends to be used as 
though it were the intrinsic quantity, rather than the one 
that has been broadened by measurement error. Our fit to 
the intrinsic distribution, based on cmodel magnitudes (from 
Table lB2l for log 10 M t /M e > 10.5) gives stellar mass den- 
sities in objects more massive than (1,2,3) x 10 Mq of 
(8.84,3.00,1.03) x lO 7 M Mpc~ 3 . This corresponds to an 
extra ~ (15,40, 109) percent more than from the Bell et al. 
fit. 

Another way to express this difference is in terms of 
the mass scale at which the integrated comoving number 
density of objects is 10 -5 Mpc -3 . For our fits to the in- 
trinsic cmodel based counts, this scale is 3.98 x 1O 11 M0 
(4.27 x 10 n Af Q for the observed fit), whereas for Bell et 
al. (2003) it is 3.31 x 1O U M . For 10" 5 5 Mpc" 3 , these 
scales are 5.25 x lO n M (5.75 x lO n M for the observed 
fit) and 4.07 x lO n M , respectively. Figure [23l illustrates 
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Figure 22. Stellar mass functions estimated from cmodel magnitudes. Top: Filled black circles show M* estimated from g — r colors 
(equation [g). Open magenta and green triangles show Af*Model an£ i Af*LRG from Blanton & Roweis (2007). Solid and dashed red 
lines show our fits to equations and lllOU with parameters reported in Table [4] for log 10 M*/Mq > 10.5 estimated from cmodel 
magnitudes. For comparison solid cyan line shows our fit to equation {9} with parameters reported in Table [4] for log 10 M*/Mq > 10.5 
estimated from Petrosian magnitudes (as in Figure [20l . Solid grey curve shows the fit reported by Bell et al. (2003) (transformed to 
Hq = 70 km s~ 1 Mpc _1 ). Bottom: Same as top panel, but now all quantities have been normalized by the Bell et al. (2003) fit. 
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Table 4. Top two rows: Parameters of <j>(L r ) (fit to M r < -17.5) and <f>{M„) (fit to log 10 M*/M > 8.5) derived from fitting 
equations (0 (in brackets) and dlOI I to the observed counts based on Petrosian magnitudes. Bottom two rows: Parameters of <j>(L r ) (fit 
to M r < —20) and (j>[M») (fit to log 10 M*/Mq > 10.5). These second set of fits better reproduce the high luminosity and mass end. 



Sample </>,/10 _2 Mpc- 3 X„ a /3 p x 



All L/1O 9 L (8.427) 8.749 ± 4.228 (15.77) 16.04 ± 2.18 (0.08) 0.08 ± 0.05 (0.827) 0.833 ± 0.036 0.128 
All M„/1O 9 M (5.620) 5.886 ± 1.839 (23.20) 25.86 ± 5.85 (0.14) 0.13 ± 0.05 (0.616) 0.654 ± 0.034 0.289 



All L/1O 9 L (1.693) 1.707 ± 0.432 (7.32) 7.56 ± 2.05 (0.55) 0.54 ± 0.16 (0.698) 0.705 ± 0.040 0.117 
A11M*/1O 9 M (0.888) 0.857 ± 0.097 (0.94) 0.95 ± 0.35 (1.39) 1.50 ± 0.16 (0.410) 0.421 ± 0.016 0.247 
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Figure 23. Comparison of various estimates of the cumulative 
number density of objects more massive than M, (top), and the 
stellar mass density in these objects (bottom). All measurements 
were transformed to Ho = 70 km s Mpc - 1 and Chabrier (2003) 
IMF. 



these differences graphically. Note that the stellar masses 
used by Panter et al. and Li & White were obtained using 
Petrosian magnitudes. To account for the difference between 
Petrosian and model luminosities one could add ~ 0.05 dex 
to their values of M* (strictly speaking, to those objects with 
C r > 2.86) - we have not applied such a shift. 
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Figure 24. Cumulative number density for comparison with 
Fig. 3 in Bczanson et al. (2009). Most of the curves are as labeled 
in top panel of Figure [23] The solid magenta line shows the cumu- 
lative number density from Marchesini et al. (2009) for galaxies 
at 2 < z < 3. The magenta dashed lines shows the same result 
but assuming that the quiescent fraction of galaxies at z ~ 2.5 is 
0.5. Dashed black line shows the number density of galaxies with 
Af» > 10 u Mq from the above fit as in Bezanson et al. (2009). 
Solid black line shows the result of having equal mass mergers. 



5.3 On major dry mergers at high masses 

The increase in z ~ counts at high masses matters greatly 
in studies which seek to constrain the growth histories of 
massive galaxies by comparing with counts at z ~ 2. Fig- 
ure [24] shows a closer-up view of the top panel of Figure [23] 
The plot is in the same format as Figure 3 in Bezanson et 
al. (2009). We have added a (solid magenta) curve showing 
the cumulative counts at 2 < z < 3 from Marchesini et al. 
(2009), shifted to our Chabrier IMF by subtracting 0.05 dex 
from their M* values. The dashed curve below it shows the 
same counts shifted downwards by a factor of two, to reflect 
the fact that only perhaps half of the galaxies at z ~ 2.5 are 
quiescent. 

Bezanson et al. (2009) argue that models in which the 
high redshift objects change their sizes but not their masses 
by the present time (e.g. Fan et al. 2008) lie well below the 
z = counts. Because this results in an order of magnitude 
fewer counts than observed at z = 0, such models, while 
viable, do not represent the primary growth mechanism of 
massive galaxies. Other models invoke minor (dry) merg- 
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Table 5. Top two rows: Parameters of 0(M») fit to log 10 M„/Mq > 8.5 (top) and > 10.5 (bottom) derived from fitting equations (in 
brackets) and dlOI I to the observed counts based on cmodel magnitudes and a Salpeter IMF for elliptical galaxies, an offset of -0.05 from 
the Salpeter IMF for SOs and an offset of -0.25 for the remaining galaxies. Bottom two rows: Similar to top two rows but for 0(Md yn ), 
where M dyn = 5R e o- 2 /G. 



Sample 


0»/lO- 2 Mpc- 3 


M„/10 9 Mq 


a 


/3 


p„/10 9 M Mpc" 3 


A IMF 


(35.196) 130.824 ± 76.761 


(69.09) 75.52 ± 7.01 


(0.01) 0.01 ± 0.02 


(0.657) 0.700 ± 0.024 


0.382 


A IMF 


(1.797) 1.958 ± 0.907 


(38.99) 49.18 ± 15.89 


(0.27) 0.23 ±0.12 


(0.595) 0.647 ± 0.046 


0.365 


M dyn 


(6.066) 6.194 ±2.699 


(17.85) 21.10 ±9.71 


(0.19) 0.18 ±0.10 


(0.485) 0.512 ± 0.041 


0.617 




(2.135) 1.757 ±0.474 


(1.28) 0.45 ± 0.28 


(0.82) 1.12 ±0.22 


(0.361) 0.337 ±0.018 


0.581 



ers (e.g., Bernardi 2009). If every one of the objects with 
M > 1O 11 M0 at z ~ 2 merged with other objects of much 
smaller mass, then the abundance of these objects would not 
change, but their masses would: the expected evolution of 
the population with M > 10 11 Mq at z ~ 2 is shown by the 
horizontal shaded region. The fractional mass increase by a 
minor merger is expected to lead to a size increase that is 
larger by a factor of two (e.g. Bernardi 2009). The observed 
size change suggests that the masses have not increased by 
more than a factor of about two: this is the vertical dashed 
line labeled 'minor mergers'. The horizontal shaded region 
intersects this vertical line at abundances which are about 
a factor of five smaller than our z = counts, so this model 
is also viable. On the other hand, if every one of the ob- 
jects with M > 1O 11 M0 at z ~ 2 merged with another 
of the same mass - a major dry merger - then this would 
shift the counts downwards and to the right, as shown by 
the shaded curved region. In this case, the fractional mass 
and size changes are equal, so the observed size increase re- 
quires mass growth by a factor of five. This is the vertical 
line labeled 'equal mass mergers'. The intersection of the 
curved shaded region with this dashed line lies above previ- 
ous estimates of the z = abundances; this lead Bezanson 
et al. (2009) to conclude that major mergers could not be 
the dominant evolution mode at the massive end. While we 
believe this an overly simplistic model, here we are simply 
pointing out that our higher abundances suggest that their 
conclusion should be revisited. 



5.4 Morphological dependence of the IMF 

Recently, Calura et al. (2009) have argued that a number of 
observations are better reproduced if one assumes a Salpeter 
(1955) IMF for ellipticals and a Scalo (1986) IMF for spirals. 
(A Salpeter IMF for ellipticals is also prefered by Treu et al. 
2010, which appeared while our paper was being refereed.) 
Whereas the M*/L-color relation for a Scalo IMF is similar 
to that for the Chabrier IMF which we have been using, 
the relation for the Salpeter IMF is offset by 0.25 dex (see 
Table [2}. Since we have found a method to separate Es, SOs 
and Spirals we can incorporate such a dependence easily. 

For ellipticals, i.e., objects selected following Hyde & 
Bernardi (2009), we compute M* by adding 0.25 to the 
right hand side of equation (|6}. For SOs, i.e., objects with 
C r > 2.86 that were not identified as ellipticals, we compute 
M» by adding 0.2 to equation (|6]), since SOs are closer to 
ellipticals than to spirals. For all other objects, we use equa- 



tion (0 as before. The open red triangles in Figure [231 show 
the stellar mass function which results. Smooth curves show 
the observed fits to equation Q; the best- fit parameters are 
reported in Table [5] 

Summing up the observed counts to estimate the stel- 
lar mass density in the range logio M*/Mq > 8.6 yields 
3.92 x 1O 8 M Mpc" 3 . Integrating the intrinsic fit (see Ta- 
ble [5]) over the entire range of masses gives a similar result: 
3.82 x 10 s M Mpc" 3 . The intrinsic fit to logi O M*/M > 
10.5 gives (2.92,1.81,0.97,0.50) x 10 s M Mpc" 3 for ob- 
jects with logioM./M© above (10.5,11.0,11.3,11.5). This 
means that ~ (75, 46, 25, 13%) of the mass is in systems 
with M» > (0.3,1,2,3) x lO n M . These estimated val- 
ues of the stellar mass density (i.e. from the intrinsic fit to 
logioM./Me > 10.5), for objects with logi O M«/M above 
(11.0,11.3,11.5), are ~ (105,224,388) percent larger than 
those infered from stellar masses computed using the cmodel 
magnitudes but using the Chabrier IMF for all types (solid 
black circles in Figure (25}- 

Finally, we compare 0(M«) with (p(Mdyn) (open blue 
squares in Figure [25)) . where Md yn = 5R e a 2 /G is the dy- 
namical mass. Our fits to 0(Md yn ) are reported in Table [S] 
At low a and R the velocity dispersions and sizes are noisy. 
This makes the mass estimate noisy below ~ 4 x 1O 9 M , 
so we only show results above this mass. Using different 
IMFs for galaxies of different morphological type reduces 
the difference between the estimated value of the stellar and 
dynamical mass especially at larger masses. 

We also find this estimate of the stellar mass function 
to be in reasonably good agreement with the one computed 
by Shankar et al. (2006) based on dynamical mass-to-light 
ratios calibrated following Salucci & Persic (1999), Cirasuolo 
et al. (2005) and references therein, lending further support 
to the possibility of an Hubble-type dependent IMF. 



5.5 The match with the integrated star formation 
rate 

It has been argued that a direct integration of the cosmolog- 
ical star formation rate (SFR) overpredicts the local stellar 
mass density (see, e.g., Wilkins et al. 2008, and references 
therein) . This has led several authors to invoke some correc- 
tions, such as a time- variable IMF. We now readdress this 
interesting issue by comparing our value for p* with that 
from integrating the SFR. 
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Figure 25. Top: Solid black circles show the stellar mass function 0(M«) estimated from cmodel magnitudes and g — r colors (equation^. 
Open red triangles show 0(M*) if we use a Salpeter IMF for elliptical galaxies, the Salpeter IMF offset by —0.05 dex in M*/L for a 
given color for SOs, and the Salpeter IMF with an offset of —0.25 for the remaining galaxies (see text for details). Open blue squares 
show o(.M ( iy U ) where M<j yn = 5R e o 2 /G. Solid and dashed lines show our observed fits (equation [U computed using log 10 M*/Mq > 10.5 
and log 10 M*/Mq > 8.6, respectively. These fits are reported in Table [5] Dot-dashed green line shows the dynamical mass function of 
Shankar et al. (2006) based on dynamical mass-to-light ratios. Bottom: Same as top panel, but now all quantities have been normalized 
by the observed fit to the black circles (from Table lB2l computed for log 10 M*/Mq > 8.6). 
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Figure 26. Cosmological SFR as given by Fardal et al. (2007; 
long-dashed line), and by Hopkins & Beacom (2006; dotted line), 
compared to the dust-corrected UV data by Bouwens et al. (2009; 
solid squares). The solid line is our "corrected" Fardal et al. SFR 
fit tuned to match the data at z > 4. 



The stellar mass density at redshift z is given by 
f z rlt' 

= J dz 'j^P^ 1 ~ /"[*(*) - * 1) • ( 16 ) 

where p* is the cosmological SFR in units of 
Mq Mpc~ 3 yr -1 , and frit) is the fraction of stellar 
mass that has been returned to the interstellar medium. 
For our IMF, 



frit) = 0.05 In 1 + 



t 



3 x 10 5 yr J ' (17) 

where t is in years (Conroy & Wechsler 2009). Note that 
this results in smaller remaining mass fractions (~ 50% at 
2 ~ instead of the usual 63-70%), than assumed in most 
previous work. This will be important below. 

We specify the SFR as follows. Bouwens et al. (2009) 
have recently calibrated the SFR over the range 2 < z < 6 
using deep optical and infrared data from ACS/NICMOS 
in the GOODS fields, UBVi droput Lyman Break Galax- 
ies, and ULIRGs data from Caputi et al. (2008). The solid 
squares in Figure [551 show their measurements, decreased 
by —0.25 dex to correct from their assumed Salpeter to a 
Chabrier IMF. These values of the SFR are much lower 
than simple extrapolations of the SFRs by Hopkins & Bea- 
com (2006) and Fardal et al. (2007), shown by dotted and 
dashed lines, respectively. We also note that the Fardal et 
al. fit matches well the updated SFR recent Bouwens et al. 
estimates in the range 2.5 < z < 4. In detail, the curves 
show the parameterization of Cole et al. (2001), 

(a + bz)y 



P*iz) 



(18) 



1 + iz/c) d ' 

Fardal et al. (2007) set a = 0.0103, b = 0.088, c = 2.4, d = 
2.8,7 = 1) an d we then multiply the total by 0.708 
(—0.1 dex) to correct from the assumed diet-Salpeter to our 
Chabrier IMF. Hopkins & Beacom (2006) set a = 0.014, b = 
0.11, c = 1.4, d — 2.2,7 = 0.7 (all parameters defined for 
h = 0.7). We convert from their IMF (from Baldry & Glaze- 
brook 2003) to the Chabrier IMF by multiplying by 1.135 
(i.e., 0.055 dex, see Table[2]). Based on detailed spectral mod- 
elling, Bouwens et al. (2009) concluded that the discrepancy 



Figure 27. Comparison of the expected stellar mass density 
based on the SFR, with the measured values at a range of red- 
shifts. Filled triangle shows our determination of the intrinsic 
local stellar mass density; filled circle shows the observed value. 
Recent determinations at higher redshifts by other groups (as 
labelled) are also shown. Dotted, dashed and solid curves show 
the result of inserting the SFRs shown in the previous figure in 
equation 116> . 



is due to dust extinction for star forming galaxies in this 
redshift range being smaller then previously assumed. (But 
we note that GRB-based estimates from, e.g., Kistler et al. 
2009, suggest this is not a closed issue.) To improve the 
match with Bouwens et al., we use the Fardal et al. values 
at z < 3.65, but set a = 0.0134, b = 0.0908, c = 3.1, d = 
6.5,7 — 0.7 at z > 3.65. This is shown by the solid curve. 

The dotted, dashed and solid curves in Figure [27] show 
the result of inserting these three models for the SFR (Hop- 
kins & Beacom, Fardal et al., and Fardal et al. corrected) 
into equation (|16[1 . The gray band bracketing the solid curve 
shows the typical ~ 15% la uncertainty (estimated by 
Fardal et al. 2007) associated with the SFR fit. The Fig- 
ure also shows a compilation of estimates of the stellar mass 
density over a range of redshifts. Our own estimate of the 
local p* value is shown by the filled triangle (filled circle 
shows the observed rather than intrinsic value) ; it is in good 
agreement with the one from Bell et al. (2003) (corrected 
to a Chabrier IMF), and is slightly larger than those from 
Panter et al. (2007), and Li & White (2009). Comparison of 
the measured p*{z) values with our new estimate of the in- 
tegrated SFR shows that the measurements lie only slightly 
below the integrated SFR at all epochs, with the discrepancy 
smallest at z > 2 and at z ~ 0. Also note that K or NIR 
selected high- 2 galaxies might be missing a significant popu- 
lation of highly obscured, dust enshrouded, forming galaxies 
in the range 0.5 < z < 2. 

The improvement with respect to previous works is 
due to the combined effects of a larger recycling factor and 
smaller high- 2 SFR, both of which act to reduce the value 
of the integral (also see discussion in Shankar et al. 2006). 
Despite the good agreement with the Fardal et al. estimate, 
we note that other SFR fits (e.g., Hopkins & Beacom) yield 
substantially higher values for the local stellar mass den- 
sity. Clearly, systematic differences such as this one must be 
resolved before this issue is completely settled. 
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Figure 28. Correlation between age and size at fixed velocity dispersion (as labeled) for types E, SO, Sa, Sb, Sc-Sd, and All in the 
Fukugita et al. (2007) sample. 



6 THE AGE-SIZE RELATION 



The previous sections studied how the distribution of L, 
M*, a and R depend on morphology or concentration. The 
present section shows one example of a correlation between 
observables which is particularly sensitive to morphology. As 



a result, how one chooses to select 'red' sequence galaxies 
matters greatly. 

The age-size relation has been the subject of recent in- 
terest, in particular because, for early-type galaxies, the cor- 
relation between galaxy luminosity and size does not depend 
on age (Shankar & Bernardi 2009). This is somewhat surpris- 
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Figure 29. Correlation between age and size at fixed dynamical mass (as labeled) for types E, SO, Sa, Sb, Sc-Sd, and All in the Fukugita 
et al. (2007) sample. 



ing, because it has been known for some time that early-type 
galaxies with large velocity dispersions tend to be older (e.g., 
Trager et al. 2000; Cattaneo & Bernardi 2003; Bernardi et 
al. 2005; Thomas et al. 2005; Jimenez et al. 2007; Shankar 
et al. 2008), and the virial theorem implies that velocity 
dispersion and size are correlated. 



To remove the effects of this correlation, Shankar et al. 
(2009c) and van der Wei et al. (2009) studied the age-size 
correlation at fixed velocity dispersion. They find that this 
relation is almost flat (with a zero-point that depends on 
the velocity dispersion, of course). At fixed dynamical mass, 
however, Shankar et al. still find no relation, whereas van der 
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Wei et al. find a significant anti-correlation: smaller galaxies 
are older. What should be made of this discrepancy? 

Although both groups claim to be studying early-type 
galaxies, the details of how they selected their samples are 
different: Shankar et al. follow Hyde & Bernardi (2009); the 
results from the previous sections suggest that this sample 
should be dominated by ellipticals, van der Wei et al. use 
the sample of Graves et al. (2009): this has 0.04 < z < 0.08, 
15 < m 7i p c t < 18, i-band concentration > 2.5; and like- 
lihood of deV profile > 1.03 that of the exponential (this 
likelihood is output by the SDSS photometric pipeline), no 
detected emission lines (EW Ha < 0.3Aand Oil < 1.7); and 
spectra of sufficient S/N that velocity dispersions were mea- 
sured (following Bernardi et al. 2003). Because the require- 
ments on the profile shape are significantly less stringent 
than those of Hyde & Bernardi, one might expect this sam- 
ple to include more SOs and Sas. Moreover, recall that age 
is a strong function of morphological type (e.g. Figure O, 
so, to see if this matters, we have studied how the age-size 
relation depends on morphological type. 

Figure [28] shows this relation, at fixed cr, for the full 
Fukugita sample (bottom right) and for the different mor- 
phological types (other panels). The relation for the full 
sample is approximately flat, except at a < 100 km s . 
However, when divided by morphological type, age (at fixed 
cr) is slightly correlated with size for ellipticals, but the rela- 
tion is more flat for SOs and Sas. (The mean age is a slightly 
increasing function of a for the later types Sb-d.) Thus, the 
flatness of the age-size relation at fixed a in the full sample 
hides the fact that the relation actually depends on mor- 
phology. 

Although this is a subtle effect for the relation at fixed 
a, the dependence on morphology is much more pronounced 
when studying galaxies at fixed Md yn = 5R e a 2 /G. Fig- 
ure [29] shows that, in the full sample, this relation is flat for 
Mdyn > 1O 11,5 M0, but decreases strongly for lower masses, 
except at Md yn < 10 10 ' 5 M©, where it appears to be curved. 
The stronger dependence here is easily understood from the 
previous figure: since dynamical mass is oc Ra 2 , as one 
moves along lines of constant mass in the direction of in- 
creasing R, one is moving in the direction of decreasing a. 
In Figure [28] this means that one must step downwards by 
one bin in logo- for every 0.4 dex to the right in log_R. For 
ellipticals, age is an increasing function of size at fixed a; 
hence, the net effect of moving up and to the right (at fixed 
cr), and then stepping down to lower cr, produces an ap- 
proximately flat age-size relation for fixed Md yn - For Sas, on 
the other hand, keeping Md yn fixed corresponds to shifting 
down and to the right (at fixed a), and then stepping down- 
wards to the lower a bin; the net result is that age decreases 
strongly as size increases. What is remarkable is that the 
ellipticals show precisely the scaling with Md yn reported by 
Shankar et al. (2009c), whereas the SOs and Sa's show that 
reported by van der Wei et al. (2009). 



7 DISCUSSION 

We compared samples selected using simple selection algo- 
rithms based on available photometric and spectroscopic in- 
formation with those based on morphological information. 
Requiring concentration indices CV > 2.6 selects a mix in 



which E+SO+Sa's account for about two-thirds of the ob- 
jects; if C r > 2.86 instead, then two-thirds of the sample 
comes from E+S0s; whereas Es alone account for more than 
two-thirds of a sample selected following Hyde & Bernardi 
(2009) (Figures QT] and [H and Table©. E's alone account 
for about 40%, 50% and 75% of the total stellar mass in 
samples selected in these three ways. 

The reddest objects at intermediate luminosities or stel- 
lar masses are edge-on disks (Figure fC?)) . As a result, samples 
selected on the basis of color alone, or cuts which run parallel 
to the red sequence are badly contaminated by such objects. 
However, simply adding the additional requirement that the 
axis ratio b/a > 0.6 is an easy way to remove such red edge- 
on disks from the 'red' sequence; the resulting sample is sim- 
ilar to requiring CV > 2.86. This may provide a simple way 
to select relatively clean early-type samples in higher red- 
shift datasets (e.g. DEEP2, zCosmos). Our measurements 
provide the low redshift benchmarks against which such fu- 
ture higher redshift measurements can be compared. 

We showed how the distribution of luminosity, stellar 
mass, size and velocity dispersion in the local universe is 
partitioned up amongst different morphological types, and 
we compared these distributions with those based on sim- 
ple selection algorithms based on available photometric and 
spectroscopic information (Figures ll4H17|) . We described our 
measurements by assuming that the intrinsic distributions 
have the form given by equation @. We showed how mea- 
surement errors bias the fitted parameters (equation I10[) . 
and used this to devise a simple method which removes this 
bias. The results, which are reported in tabular form in Ap- 
pendix^ show that ellipticals contain ~20% of the luminos- 
ity density and 25% of the stellar mass density in the local 
universe, and have mean sizes of order 3.2 kpc. Including 
SOs increases these numbers to 33% and 40%; adding Sas 
results in further increases, to 50% and 60% respectively. 
These numbers are in broad agreement with those from the 
Millennium Galaxy Survey of about 10 4 objects in 37.5 deg 2 . 
Driver et al. (2007) report that 15 ± 5% of the stellar mass 
density is in ellipticals, and adding bulges increases this to 
44 ±9%. 

Our stellar mass function has more massive objects than 
other recent determinations (e.g. Cole et al. 2001; Bell et 
al. 2003; Panter et al. 2007; Li & White 2009), similarly 
shifted to a Chabrier (2003) IMF (Figures [20l and |22|) . The 
mass scale on which the discrepancy arises is of order where 
some previous work had only a handful of objects - our 
substantially larger volume is necessary to provide a more 
reliable estimate of these abundances. Using stellar masses 
estimated from cmodel luminosities, which are more reliable 
than Petrosian luminosities at the large masses where the 
discrepancy in <j>(M*) is largest, gives stellar mass densities 
in objects more massive than (1,2,3) x 10 Mq that are 
larger by more than ~ (20, 50, 100) percent compared to 
Bell et al. (2003) (Figure [25) . 

This analysis required that we study the sytematic dif- 
ferences between the stellar mass estimates based on g — r 
color (our equation [5] following Bell et al. 2003), colors in 
multiple bands (Blanton & Roweis 2007), and on spectral 
features (Gallazzi et al. 2005). (See Gallazzi & Bell 2009, 
which appeared while our work was being refereed, for a 
discussion of the pros and cons of these various approaches, 
and of the accuracy to which stellar masses can currently be 
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derived.) The g — r and Gallazzi et al. estimates are gener- 
ally in good agreement (Figure I Alp , although the spectral 
based estimates suffer slightly from aperture effects which 
are complicated by the magnitude limit of the survey (Fig- 
ures [X2] \KE\ and [20} . The Blanton et al. estimates are in 
good agreement with the other two provided one uses LRG- 
based templates to estimate masses at the most massive end 
(Figures IA31 IA4 l and I22 [) . At lower masses, some combina- 
tion of the LRG and other templates is required. Ignoring 
the LRG templates altogether (e.g. Li & White 2009) re- 
sults in systematic underestimates of as much as 0.1 dex 
or more (Figure IA3|) . severely compromising estimates of 
the number of stars currently locked up in massive galaxies 
(Figure [20]l . If we compare our estimate of the stellar mass 
density in objects more massive than (1, 2, 3) x 10 11 M© with 
those from the Li & White (2009) fit, then our values are 
~ (140,230,400) percent larger. 

Allowing more high mass objects means that major dry 
mergers may remain a viable formation mechanism at the 
high mass end (Figure |2"4"|). It also relieves the tension be- 
tween estimates of the evolution of the most massive galaxies 
which are based on clustering (which predict some merg- 
ing, and so some increase in stellar mass; Wake et al. 2008; 
Brown et al. 2008) and those based on abundances (for which 
comparison of high redshift measurements with the previous 
z ~ measurements indicated little evolution; Wake et al. 
2006; Brown et al. 2007; Cool et al. 2008). This discrepancy 
may be related to the origin of intercluster light (e.g. Skibba 
et al. 2007; Bernardi 2009); our measurement of a larger lo- 
cal abundance in galaxies reduces the amount of stellar mass 
that must be stored in the ICL. 

It has been argued that a number of observations are 
better reproduced if one assumes a different IMF for ellipti- 
cal and spiral galaxies (e.g. Calura et al. 2009). We showed 
that this acts to further increase the abundance of mas- 
sive galaxies (Figure [23]). and reduces the difference between 
stellar and dynamical mass, especially at larger masses. At 
M, > lO n M , the increase due to the change in IMF is a 
factor of two with respect to models which assume a fixed 
IMF. 

If we sum up the observed counts to estimate the stellar 
mass density in the range 8.6 < log 10 M* < 12.2 Mq (M* 
from equation [6] using cmodel magnitudes), then the result 
is 3.05 x 10 s Mq Mpc -3 . Using our fit to the observed distri- 
bution (values between round brackets in Table |B2} gives a 
similar value (3.06 x 10 s Mq Mpc -3 ) and a slightly smaller 
value if one uses the intrinsic fit (2.89 x 10 8 Mq Mpc -3 , see 
Table |B2} . Our values are ~ 15% and 30% larger than those 
reported by Panter et al. (2007) and Li & White (2009), 
respectively. If we allow a type dependent IMF, the total 
stellar mass density increases by a further 30%. 

However, although our stellar mass function has more 
M* > 10 11 Mq objects than other recent determinations, 
our estimate of the total stellar mass density is similar to 
that measured by Bell et al. (2003). It is about 20% smaller 
than the value reported by Driver et al. (2007) (once shifted 
to the same IMF, for which we have chosen Chabrier 2003; 
see Table [2]). This is because differences at the mid/faint 
end contribute more to the total stellar density than the 
difference we measured at the massive end. 

It has been suggested that direct integration of the cos- 
mological star formation rate overpredicts the total local 



estimate of the stellar mass density (see, e.g., Wilkins et al. 
2008, and references therein). However, we showed that re- 
cent determinations of the recycling factor (equation ll6|) and 
the high- 2 star formation rate (Figure I26|) result in better 
agreement (Figure I27[) . This is because the former yields 
smaller remaining masses, and the latter produces fewer 
stars formed in the first place. 

Our measurements also show that the most luminous 
or most massive galaxies, which one might identify with 
BCGs, are less concentrated and have smaller b/a ratios, 
than slightly less luminous or massive objects (Figures [3] 
[5] and Hip . Their light profile is also not well represented 
by a pure deVaucoleur law. This is consistent with results 
in Bernardi et al. (2008) and Bernardi (2009) who suggest 
that these are signatures of formation histories with recent 
radial mergers. In this context, note that we showed how to 
define seeing-corrected sizes, using quantities output by the 
SDSS pipeline, that closely approximate deVaucouleur bulge 
+ Exponential disk decompositions (equations[5]and|3J . Our 
cmodel sizes represent a substantial improvement over Pet- 
rosian sizes (which are not seeing corrected) and pure deV 
or Exp sizes (Figures [TJ and [3}. 

And finally, our study of the age-size correlation re- 
solves a discrepancy in the literature: whereas Shankar et 
al. (2009c) report no correlation at fixed Md yn , van der Wei 
et al. (2009) report that larger galaxies tend to be younger. 
We showed that ellipticals follow the scaling reported by 
Shankar et al. scaling, whereas SOs and Sas follow that of 
van der Wei et al. (Figure |2"9"} . suggesting that Shankar et al. 
select a sample dominated by galaxies with elliptical mor- 
phologies, whereas van der Wei et al. include more SOs and 
Sas. These conclusions about the differences between the 
samples are consistent with how the samples were actually 
selected. 

Since van der Wei et al. use their measurements to con- 
strain a model for early-type galaxy formation, this is an 
instance in which having morphological information mat- 
ters greatly for the physical interpretation of the data. Our 
results indicate that models of early-type galaxy formation 
should distinguish between ellipticals and SOs because, in 
the projection of the age-size-mass correlation shown in Fig- 
ure !29l the SOs and Sas are very different from the other mor- 
phological types. Whether the smaller sizes for older SOs are 
due to the gradual stripping away of a younger disk is an 
open question. 

van der Wei et al. (2009) use their observation that the 
age-size relation at fixed a is flat to motivate a model in 
which early-type galaxy formation requires a critical veloc- 
ity dispersion (which they allow to be redshift dependent). 
The same logic applied here suggests that while this may be 
reasonable for SOs or Sas, it is not well-motivated for ellipti- 
cals (Figure [28}. However, it might be interesting to explore 
a model in which elliptical formation requires a critical (pos- 
sibly redshift dependent) dynamical mass rather than ve- 
locity dispersion (which may be redshift dependent). This is 
interesting because, in hierarchical models, the phenomenon 
known as down-sizing (e.g., Cowie et al. 1996; Heavens et al. 
2004; Sheth et al. 2006; Jimenez et al. 2007) is then easily 
understood (Sheth 2003). 
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APPENDIX A: COMPARISON OF STELLAR 
MASS ESTIMATES 

This section shows results from a study of how the different 
M» estimates compare with one another. We first show re- 
sults for galaxies that have CV > 2.86 for two reasons. First, 
because this is expected to provide a sample dominated by 
massive galaxies, and these are expected to be more ho- 
mogeneous. And second, we are particularly interested in 
studying the high mass end of the stellar mass function, be- 
cause this is the end which is first observed in higher redshift 
datasets. 

Figure IA1I shows that the M*Beii and M*Gaiiazzi esti- 
mates are in good agreement, although M*Bcii tends to be 
larger at M.bcII > 10 M©. The sharp downturn at small 
masses is due to aperture effects. This is because M*b c ii is 
based on the model color; for otherwise similar galaxies, this 
samples the same fraction of a galaxy's light (that within R e ) 
whatever its distance. On the other hand, M*Gaiiazzi is based 
on the light which enters the SDSS fiber (radius 1.5 arcsec); 
this samples a distance dependent fraction of the light of a 
galaxy. For small low mass galaxies which have significant 
disks, the fiber samples the bulge of the nearby objects, and 
the disks of the more distant ones. Since M,/L is larger for 
bulges than disks, the more distant objects of a given M*b c ii 



will have smaller M*Gaiiazzi- Here, of course, we have explic- 
itly selected against disky galaxies, but color gradients will 
produce the same effect. 

Figure [A2l shows this effect explicitly: at fixed M*Beii, 
the ratio M t Gaiiazzi/Af*Boii decreases with increasing dis- 
tance (this is a more dramatic effect for the lower mass 
systems in samples which include disks - see Figure IA5I) . 
The panel on the right, which shows the angular size for 
objects of a fixed M^Beii increasing at high z, illustrates a 
curious selection effect (the upturn at higher redshift for the 
middle stellar mass bin, which is not present if when we 
study bins in luminosity). For a given M.bcII, the highest 
redshift objects in a magnitude limited survey will be bluer; 
since color is primarily determined by velocity dispersion, 
these objects will have smaller than average velocity disper- 
sions, and larger sizes. Although this is a small effect for 
the massive objects we have selected here, it is more dra- 
matic in samples which include disks (as we show shortly; 
see Figure [A"5"t . 

Figure IA3I shows how these estimates compare with 
those from Blanton & Roweis (2007). The top left panel 
shows that M^Modei is ~ 0.04 dex larger than Af*p otro (top 
left panel) . Because the two estimates are based on the same 
mass-to-light ratio, this offset is entirely due to the fact 
that Petrosian magnitudes are 0.1 mags fainter than (the 
more realistic) model magnitudes. The middle left panel 
shows that M.Beii/Aif.Modci is approximately independent 
of mass for g — r, and strongly mass dependent for r — i. 
But perhaps most importantly, note that M.Modci is typ- 
ically 0.1 dex smaller than the g — r based M«Beli' This 
offset would have been even larger if we had compared to 
M*p e tro - the fairer comparison, because our M*Beii here 
is based on the Petrosian magnitude. We argue shortly 
that this offset reflects the fact that the templates used 
to estimate M*Modei are not appropriate for large masses. 
In contrast, M»G a ii azz i/M*Modei is strongly mass dependent 
at small masses - this is because aperture effects matter 
for MtGaiiazzi but not for M*Model- The behaviour at larger 
masses is a combination of aperture and template effects, 
but note again that M*Model is about 0.1 dex low. 

The panels on the right show that A/*lrg fares much 
better. The top right panel shows that it is always substan- 
tially larger than M*Modei, and the two panels below it show 
that M.lrg is much more like M*Bell and M*Gaiiazzi than 
is M„ Model- Comparison with the corresponding panels on 
the left shows that much of the offset is removed, strongly 
suggesting that neither M*p otr o nor M*Modei are reliable at 
the high mass end. This is not surprising - massive galaxies 
are expected to be old, so a 10 Gyr template should provide 
a more accurate mass estimate than ones that are restricted 
to ages of ~ 7 Gyrs or less. This will be important in Sec- 
tion [5] At the low mass end, M*lrg is larger than either 
of the other estimates, suggesting that it overestimates the 
true mass. However, at these smaller masses, using M*Modei 
exclusively produces masses which lie below M*Gaiiazzi and 
MtBeii- This suggests that the LRG template remains the 
better description for a non-negligible fraction of low mass 
objects. Therefore, in practice, to use the Blanton & Roweis 
estimates, one must devise a method for choosing between 
M» Model and Af.LRG- 

Figure [A"4l compares the different estimates of the stel- 
lar mass for the full sample. This gives similar results to 
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Figure Al. Comparison of stellar masses computed following Bell et al. (2003) (our equations [6] (solid circles) and [7] (open squares) 
with Petrosian r-band luminosity), and Gallazzi et al. (2005), for objects with CV > 2.86. 
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Figure A2. Aperture effects on the Gallazzi et al. (2005) stellar mass estimate. Top: Redshift dependence of M t Q a ii azz ;/M,,Bcll (left) 
and the angular half-light radius (right) for objects with CV > 2.86 in three different bins of M 4Bc ll as indicated. Bottom: Similar, but 
now for a few narrow bins in luminosity. 



when we restricted to CV > 2.86, although the offsets are 
more dramatic. In this case, the selection effect associated 
with studying fixed Af.Bcii makes the aperture effect asso- 
ciated with M*Gaiiazzi appear more complex than it really 
is: the large bump at the highest redshift associated with a 
given M*Bcii is not present when one studies objects at fixed 
luminosity (compare top and bottom panels of Figure |A"5| . 
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Figure A3. Comparison of various stellar mass estimates from Blanton &; Roweis (2007), with M^Bell and M*Gallazzi for objects with 
CV > 2.86. 
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Figure A4. Comparison of our various stellar masses in the full sample (i.e., no cut on concentration index). 





Figure A5. Same as Figure \K2\ but now for the full sample (i.e. no cut on concentration index). 
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APPENDIX B: TABLES 

This Appendix provides tables which summarize the results 
of fitting equation (|10p to the measured luminosity, stel- 
lar mass, size and velocity dispersion distributions shown 
in Figures I14H17I These measurements were obtained us- 
ing cmodel magnitudes and sizes. We described our mea- 
surements by assuming that the intrinsic distributions have 
the form given by equation ©. The values between round 
brakets show the parameters obtained fitting equation ([9]) 
to the data ignoring measurement errors. 
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Table Bl. Best-fit parameters of equation l|10[ l to the measured r— band cmodel luminosity function <j>(L). Values between round brakcts 
show the parameters obtained fitting equation J9j to the data ignoring measurement errors. 



Sample 


</>„/icr 2 Mpc- 3 


L./10 9 Lq 


a 


P 


Pl/10 9 L Mpc- 3 


HB09 


(0.095) 0.095 ±0.005 


(10.69) 10.99 ±4.06 


(1.38) 1.37 ±0.18 


(0.769) 0.776 ±0.071 


0.024 


CI > 2.86 


(0.174) 0.174 ±0.010 


(6.96) 7.21 ± 2.52 


(1.38) 1.37 ± 0.16 


(0.692) 0.698 ± 0.051 


0.038 


CI > 2.6 


(0.382) 0.382 ± 0.022 


(3.16) 3.28 ± 1.30 


(1.32) 1.31 ± 0.16 


(0.583) 0.588 ±0.039 


0.061 


All 


(4.123) 4.196 ± 1.693 


(12.21) 12.48 ±3.42 


(0.20) 0.19 ± 0.10 


(0.728) 0.734 ± 0.051 


0.135 


All (M r < -20) 


(1.227) 1.197 ± 0.240 


(2.08) 1.88 ± 0.86 


(1.07) 1.12 ± 0.24 


(0.540) 0.533 ±0.035 


0.125 


F07-E 


0.13 ±0.11 


72.10 ±34.38 


0.320 ± 0.444 


1.752 ± 0.750 


0.022 


F07-S0 


1.05 ±0.65 


64.05 ± 10.15 


0.065 ± 0.046 


1.462 ± 0.203 


0.038 


F07-Sa 


0.80 ± 0.07 


27.19 ± 1.15 


0.254 ± 0.020 


0.944 ± 0.023 


0.058 


F07-Sb 


17.22 ± 1.46 


36.01 ± 1.35 


0.014 ± 0.001 


1.079 ± 0.030 


0.084 


F07-Scd 


23.48 ± 1.93 


28.37 ± 1.04 


0.017 ± 0.001 


0.980 ± 0.023 


0.112 


F07-A11 


16.00 ±8.98 


25.67 ± 4.87 


0.027 ± 0.016 


0.934 ± 0.086 


0.116 



Table B2. Best-fit parameters of equation (1 1 I t to the measured stellar mass function <f>(M*). 



Sample 


<A*/KT 2 Mpc" 3 


M./10 9 M Q 


a 


P 


p,/10 9 M Q Mpc~ 3 


HB09 


(0.095) 0.095 ± 0.005 


(20.25) 25.24 ±8.70 


(1.29) 1.28 ± 0.14 


(0.641) 0.696 ±0.051 


0.066 


CI > 2.86 


(0.174) 0.174 ±0.009 


(13.62) 17.62 ±6.14 


(1.26) 1.24 ± 0.13 


(0.590) 0.640 ± 0.043 


0.105 


CI > 2.6 


(0.388) 0.389 ± 0.020 


(7.49) 10.07 ± 3.89 


(1.10) 1.07 ± 0.12 


(0.522) 0.563 ± 0.036 


0.166 


All 


(5.285) 5.850 ± 3.054 


(34.88) 39.06 ± 8.81 


(0.11) 0.09 ± 0.06 


(0.650) 0.694 ± 0.040 


0.289 


All (M* > 3 x 1O 1O M ) 


(0.761) 0.672 ±0.123 


(0.14) 0.02 ± 0.01 


(1.95) 2.68 ± 0.30 


(0.342) 0.308 ± 0.010 


0.261 


F07-E 


0.09 ±0.04 


158.43 ± 115.34 


0.54 ±0.54 


1.31 ± 0.53 


0.062 


F07-S0 


1.11 ±0.66 


206.53 ± 35.47 


0.05 ± 0.04 


1.44 ± 0.25 


0.104 


F07-Sa 


1.70 ±2.52 


144.73 ± 45.50 


0.07 ±0.11 


1.18 ± 0.23 


0.151 


F07-Sb 


1.46 ± 1.31 


90.24 ±46.21 


0.15 ± 0.18 


0.96 ± 0.21 


0.203 


F07-Scd 


4.09 ±5.46 


81.34 ±30.52 


0.07 ±0.11 


0.92 ± 0.15 


0.248 


F07-A11 


6.76 ±9.92 


84.20 ± 24.84 


0.04 ±0.07 


0.93 ±0.13 


0.251 



Table B3. Best-fit parameters of equation l|10| l to the measured r— band cmodel size function <j>(R e ). 



Sample 


</>„/10- 2 Mpc- 3 


R*/kpc 


a 


P 


(R)/kpc 


HB09 


(0.093) 0.093 ± 0.006 


(0.0001) 0.0003 ± 0.0002 


(8.70) 9.29 ±0.45 


(0.293) 0.294 ±0.007 


3.25 


CI > 2.86 


(0.173) 0.173 ±0.011 


(0.0001) 0.0002 ± 0.0002 


(9.82) 10.73 ± 0.53 


(0.294) 0.290 ± 0.007 


3.19 


CI > 2.6 


(0.400) 0.400 ± 0.024 


(0.0001) 0.0002 ± 0.0002 


(7.58) 8.08 ±0.38 


(0.295) 0.295 ± 0.006 


2.74 


All 


(6.009) 6.040 ± 0.766 


(0.4838) 0.5809 ± 0.1773 


(1.31) 1.27 ± 0.23 


(0.688) 0.729 ± 0.045 


1.41 


All (1.5 < R e < 20 kpc) 


(2.951) 2.858 ±0.420 


(0.5523) 0.5612 ± 0.2033 


(2.22) 2.39 ±0.44 


(0.765) 0.790 ± 0.058 


2.38 


F07-E 


0.10 ± 0.02 


0.0007 ± 0.0009 


6.44 ± 1.55 


0.352 ± 0.032 


3.13 


F07-S0 


0.19 ±0.03 


0.0008 ± 0.0008 


6.07 ± 1.05 


0.356 ± 0.025 


2.78 


F07-Sa 


0.37 ±0.06 


0.0004 ± 0.0004 


6.08 ± 1.01 


0.337 ±0.021 


2.46 


F07-Sb 


0.83 ± 0.15 


0.0123 ± 0.0104 


4.24 ±0.82 


0.444 ± 0.037 


2.30 


F07-Scd 


1.55 ±0.21 


0.0062 ± 0.0042 


5.72 ±0.76 


0.431 ± 0.026 


2.80 


F07-A11 


3.92 ±0.77 


0.0049 ± 0.0033 


4.02 ± 0.64 


0.402 ± 0.024 


1.81 
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Table B4. Best-fit parameters of equation l|10[ l to the measured velocity dispersion function 0(c). 



Sample 


j, /i n— 2iv/rT-,r>— 3 
(p* i iu ivipc 


u # /km s ^ 


Oi 


a 
P 


(c } /km s ■ 


HB09 


(0.097) 0.096 ± 0.006 


(184.08) 173.41 ± 16.85 


(2.44) 2.83 ±0.41 


(2.91) 3.10 ±0.35 


149.15 


CI > 2.86 


(0.182) 0.179 ±0.010 


(177.34) 166.50 ± 17.02 


(2.17) 2.54 ±0.41 


(2.76) 2.93 ±0.33 


139.14 


CI > 2.6 


(0.663) 0.590 ±0.088 


(190.57) 175.96 ± 16.88 


(0.80) 1.06 ± 0.34 


(2.86) 2.99 ±0.35 


92.40 


All 


(2.099) 2.099 ± 0.099 


(113.78) 113.78 ± 1.06 


(0.94) 0.94 ±0.03 


(1.85) 1.85 ±0.02 


63.70 


All (<r > 125 km s" 1 ) 


(8.133) 2.611 ± 0.161 


(176.99) 159.57 ± 1.48 


(0.11) 0.41 ± 0.02 


(2.54) 2.59 ±0.04 


44.02 


F07-E 


0.11 ± 0.04 


218.27 ± 43.06 


1.53 ± 1.21 


4.47 ± 1.90 


131.56 


F07-S0 


0.19 ± 0.05 


197.35 ± 58.97 


1.66 ± 1.25 


3.49 ± 1.75 


128.40 


F07-Sa 


0.41 ± 0.13 


193.07 ± 44.05 


1.06 ±0.77 


3.35 ± 1.29 


99.32 


F07-Sb 


0.61 ± 0.25 


158.18 ± 59.10 


1.22 ± 1.14 


2.57 ± 1.05 


93.12 


F07-Scd 


2.47 ± 3.96 


180.45 ± 33.69 


0.29 ±0.54 


2.96 ±0.85 


37.48 


F07-A11 


1.23 ±0.31 


59.04 ±31.48 


2.44 ± 1.17 


1.35 ±0.30 


86.74 



